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This  report  summarizes  the  results  of  investigations  into 
the  effects  on  explosion-generated  ground  motion  of  the  MX 
system  siting  environment  in  alluvial  valleys.  The  objective 
of  the  analyses  described  in  the  report  has  been  to  develop 
a  better,  quantitative  understanding  of  the  late-time,  low- 
frequency  ground  motion  expected  in’  such  environments  in 
order  to  provide  a  firmer  basis  for  scaling  to  untested 
geologic  environments  and  attack  conditions.  The  studies  \ 
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reported  have  focused  on  two  aspects  of  the  problem  of  low-f requency 

surface  waves  which  dominate  the  ground  motion  generated  by  blasts  in 
MX  valleys:  (1)  definition  of  site  and  yield  dependence  of  ground 
motion  characteristics  at  alluvial  valley  sites,  and  (2)  possible  en¬ 
hancement  of  ground  motion  at  sites  in  alluvial  valleys  due  to  surface- 
wave  reverberation. 

With  regard  to  first  aspect,  we  have  relied  primarily  on  the 

theory  of  surface-wave  generation  by  air-blast  loading  on  an  alluvial 
valley  to  determine  scaling  relations  as  a  function  of  both  explosion 
yield  and  site  geology.  The  amplitude  of  the  Rayleigh  wave  source 
function  has  been  found  to  be  directly  proporational  to  yield  at  low 
frequencies  and  independent  of  site  geology;  while,  at  high  frequencies, 
it  behaves  as  yield  to  the  one-third  power  and  increases  with  near¬ 
surface  shear  wave  velocity.  These  relations  have  been  applied  to 
several  specific  prototype  MX  sites;  and  it  has  been  found  that  pre¬ 
dominant.  Rayleigh  wave  frequencies  at  the  sites  depend  only  on  thick¬ 
ness  and  shear  velocity  of  the  surficial  soil  layer  and  are  independ¬ 
ent  of  yield  for  yields  greater  than  a  few  hundred  kilotons.  An  anoma¬ 
lous  arrival  observed  in  the  Misers  Bluff  II-2  experiment  (the  ''800 
millisecond  anomaly")  has  been  found  to  be  inconsistent  with  simple 
linear  superposition  of  single  burst  ground  motions  and  may  correspond 
to  a  true  multiburst  effect.  The  anomalous  ground  motion  has  been 
shown  to  match  a  theoretical  model  which  considers  Rayleigh  waves  con¬ 
verging  on  the  center  of  the  source  array  from  an  axisymmetric ,  im¬ 
pulsive  surface  source,  acting  at  a  radius-coincident  with  that  of  the 
charge  locations. 

In  the  valley  reverberation  element  of  the  study,  emphasis  has  been 
on  analysis  of  observational  data,  primarily  the  ground  motions  re¬ 
corded  during  the  Misers  Bluff  II  experiment  at  Planet  Ranch  in 
Arizona  in  1978  .  Observations  from  Misers  Bluff  II  were  found  .to  be 
consistent  with  previous  experimental  results  as  well  as  theoretical 
calculations  which  indicate  that  late-time,  low-frequency  ground 
motion  in  alluvial  valley  environments  are  modified  by  the  presence 
of  valley  boundaries.  In  particular,  it  has  been  found  that  duration 
of  ground  motion  in  such  environments  may  be  significantly  prolonged. 
Differences  between  the  late  time  ground  motions  for  the  MBII-1  and 
MB I I - 2  tests  suggest  that  surface-wave  reverberation  has  a  complex 
dependence  on  source  and  receiver  locations  within  the  valley  struc¬ 
ture  . 
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I. 


INTRODUCTION 


Both  experimental  and  theoretical  studies  that  have 
been  conducted  during  the  past  several  years  have  confirmed 
that  the  potential  ground  motion  environment  for  the  MX  sys¬ 
tem  may  be  considerably  different  than  that  for  which  pre¬ 
vious  strategic  systems  have  been  designed.  In  particular, 
the  Multiple  Protective  Structure  concept  together  with  the 
site  requirements  have  combined  to  increase  the  relative  im¬ 
portance  of  the  late-time,  longer-period  ground  motions  as¬ 
sociated  with  elastic,  surface  wave  propagation  at  the  site. 
That  is,  the  system  will  be  sited  in  a  valley  and,  under  a 
multiple  attack  scenario,  untargeted  points  can  be  expected 
to  experience  a  significant  ground  motion  environment  as  the 
result  of  the  combination  of  effects  originating  from  attacks 
on  surrounding  aimpoints.  Furthermore,  reflections  of  out¬ 
going  energy  from  the  valley  boundaries  can  be  expected  to 
complicate  the  ground  motion  environment  within  the  valley, 
particularly  at  late  times. 

The  objective  of  the  analyses  described  in  this  re¬ 
port  has  been  to  develop  a  better,  quantitative  understanding 
of  these  late-time,  long-period  ground  motions  in  order  to 
provide  a  firmer  basis  for  scaling  to  untested  geologic  en¬ 
vironments  and  attack  conditions.  Two  specific  problem  areas 
have  been  addressed  during  the  past  year.  The  first  concerns 
the  definition  of  the  site  and  yield  dependence  of  the  char¬ 
acteristics  of  the  low-frequency  ground  motion  at  a  variety 
of  prototype  MX  sites.  Progress  in  this  study  area  is  sum¬ 
marized  in  Section  II  which  includes  a  detailed  discussion  of 
scaling  relations  for  airb last-induced  Rayleigh  waves  with 
applications  to  five  specific  sites  of  potential  interest,  as 
well  as  preliminary  analysis  of  a  new,  multiburst  surface 
wave  phenomenon  observed  on  the  Misers  Bluff  II-2  experiment. 
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The  second  problem  area  treated  in  this  report  con¬ 
cerns  an  assessment  of  the  possibility  that  the  MX  ground 
motion  environment  may  be  significantly  altered,  and  in  par¬ 
ticular  enhanced,  as  a  result  of  siting  the  system  in  a  val¬ 
ley.  During  the  past  several  years  a  number  of  specific  ex¬ 
perimental  and  theoretical  investigations  have  been  conducted 
in  an  attempt  to  define  the  effects  of  valley  boundaries  on 
the  ground  motion  environment  to  be  expected  at  target  points 
inside  the  valley.  A  brief  summary  of  the  results  of  these 
earlier  studies  is  presented  in  Section  III  together  with  a 
detailed  analysis  of  the  valley  reverberation  data  recorded 
from  the  Misers  Bluff  II  experiments.  This  is  followed  in 
Section  IV  by  a  summary,  together  with  a  statement  of  con¬ 
clusions  summarizing  our  current  state  of  knowledge  with  re¬ 
gard  to  these  issues. 
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II.  RAYLEIGH  WAVES  FROM  NEAR-SURFACE  EXPLOSIONS 


2 . I  INTRODUCTION 

As  a  result  of  research  conducted  during  the  past  sev¬ 
eral  years,  much  progress  has  been  made  in  defining  the  char¬ 
acteristics  of  the  ground  roll  component  of  the  motion  pro¬ 
duced  by  near-surface  explosions.  In  particular,  it  has  been 
demonstrated  (Murphy,  1978;  Auld  and  Murphy,  1979)  that  these 
low-frequency  arrivals  are  predominantly  elastic  Rayleigh 
waves  driven  by  the  airblast  loading  acting  on  the  surface 
outside  the  region  of  strong  nonlinear  interaction  (i.e.,  the 
crater) .  Moreover,  it  has  been  shown  that  these  Rayleigh 
waves  can  be  a  very  significant  component  of  the  ground  mo¬ 
tion,  particularly  at  ranges  where  the  peak  airblast  over¬ 
pressure  is  less  than  a  few  hundred  psi  (pounds  per  square 
inch) .  For  example,  Figure  1  shows  a  comparison  between  the 
vertical  particle  velocity  waveform  measured  at  a  range  of 
355  m  from  the  Pre-Mine  Throw  100  ton  surface  HE  experiment 
and  the  theoretically  predicted  Rayleigh  wave  motion  (Murphy, 
197  8)  .  It  can  be  seen  that  the  surface  wave  component  domi¬ 
nates  the  motion  at  this  range  and,  furthermore,  that  it  can 
be  simulated  with  remarkably  good  accuracy  using  the  selected 
theoretical  model.  A  similar  capability  has  also  been  demon¬ 
strated  for  the  HE  explosions  at  other  sites  (e.g.,  Pre-Dice 
Throw,  Murphy,  1978;  Misers  Bluff  II,  Murphy  and  Bennett, 

1979)  indicating  that  the  dependence  of  the  surface  wave 
characteristics  on  site  geology  can  be  accounted  for  within 

the  framework  of  a  quantitative  model. 

/ 

The  significance  of  these  findings  with  regard  to  the 
present  investigation  lies  in  the  fact  that  the  theory  gov¬ 
erning  the  generation  and  propagation  of  elastic  Rayleigh 
waves  is  well  developed  and  thus  can  be  used  to  help  define 
scaling  laws  for  the  low-frequency  components  of  the  ground 
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wave  motion,  Pre-Mine  Throw,  vertical  component, 
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motion  environment.  In  particular,  the  results  of  studies 
conducted  to  date  suggest  that  it  should  be  possible  to  esti¬ 
mate  the  dominant  frequency  and  relative  duration  and  ampli¬ 
tude  level  of  the  low-frequency  ground  motion  to  be  expected 
at  different  MX  prototype  sites  as  a  function  of  site  geology 
and  explosion  yield.  In  the  following  sections,  the  theo¬ 
retical  basis  for  the  surface  wave  model  will  be  reviewed  and 
applied  to  a  preliminary  assessment  of  some  of  the  problems 
associated  with  the  definition  of  the  low- frequency  ground 
motion  environment  for  the  MX  system. 

2.2  DESCRIPTION  OF  THE  MODEL 

The  problem  of  interest  here  concerns  the  elastic  re¬ 
sponse  of  a  multilayered  half space  to  a  propagating,  axi- 
symmetric  normal  load  acting  on  the  surface  (i.e.,  an  approxi¬ 
mation  to  an  airblast) .  The  geometry  of  the  problem  is  il¬ 
lustrated  in  Figure  2  which  shows  a  sequence  of  constant 
thickness  (d^)  layers  overlying  a  semi- inf inite  halfspace 
being  acted  on  by  a  surface  source  centered  on  the  origin. 

As  a  first  approximation,  the  layers  are  taken  to  be  per¬ 
fectly  elastic  and  are  characterized  by  their  Lame  constants 
(\,u)  and  density  (p).  The  complete  analytic  solution  to 
this  boundary  value  problem,  representing  the  total  response 
of  the  multilayered  half space,  is  very  complicated  and  not 
yet  fully  developed.  However,  the  surface  wave  portion  of 
the  solution  is  simple  enough  that  the  exact  result  can  be 
written  down  in  compact  form  using  the  matrix  formulation  of 
the  problem  pioneered  by  Haskell  (1953)  and  Harkrider  (1964) . 
In  fact,  for  the  airblast  loading  condition  being  considered 
here,  it  can  be  shown  (Harkrider,  et  al. ,  1974;  Bache,  et  al . 
1977)  that  the  vertical  (W)  and  radial  (U)  components  of  the 
Rayleigh  wave  contribution  to  the  motion  at  horizontal  range 
r  outside  the  source  region  can  be  written  in  the  form: 
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w.  (<o,r)  =  -iirA_  (<u)  H^2)  (k.r)  S .  (a) 

J  w  j  j 


Oj  (u,r)  *  c j  (cu) 


h{2) (k.r) 

,  («Q)  4^  3 


Hj2T(k.r) 


W. (w,r) 


(2-1) 


where  subscript  j  denotes  the  mode  number,  Ar.  (oj)  is  the 

( 2 )  *  J 

Rayleigh  wave  site  response  function,  (k^r)  is  the  Hankel 

function  propagation  term  for  wavenumber  k,  (i  =  0,1),  (w) 

is  the  ellipticity  function  and  S ^  (w)  is  the  Rayleigh  wave 
source  function  given  by: 


i:(“)  ’I 


p  (r"  ,cd)  JQ  (k^r")  r'dr' 


(2-2) 


where  p(r,w)  is  the  Fourier  transform  of  the  airblast  load 
acting  over  a  surface  disk  extending  from  r^  to  r2  (°£r2_£ 
r2  £r)  and  J^k^r)  is  the  Bessel  function  kernel. 

Equations  (2-1)  and  (2-2)  define  the  solution  for  the 
Fourier  transform  of  the  displacement  at  distance  r  due  to 
Rayleigh  waves  diverging  from  a  central  source  region,  and 
thus  represents  the  complete  solution  for  far-field  observa¬ 
tion  points.  However,  in  the  current  application,  the  situ¬ 
ation  is  complicated  by  the  fact  that  the  observation  points 
are  at  near-field  distances  and  may,  in  fact,  lie  inside  the 
loaded  region.  In  such  cases,  it  is  necessary  to  also  con¬ 
sider  the  contributions  to  the  motion  arising  from  Rayleigh 
waves  converging  on  the  origin  from  the  exterior  source 
region.  The  source/receiver  geometries  for  these  diverging 
and  converging  cases  are  compared  in  Figure  3.  The  solution 
for  the  Rayleigh  wave  component  of  the  displacement  in  the 
converging  wave  case  can  be  deduced  from  an  analysis  of  the 
simpler  halfspace  problem  (Bache,  1979;  Murphy  and  3ache, 
1979) .  Consider  the  case  in  which  an  axi symmetric  impulsive 
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Figure  3 .  Comparison  of  diverging  and  converging  Rayleigh 
wave  geometry  for  axisymmetric  loads  acting  on 
the  surface  of  a  multilayered  halfspace. 
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surface  load  is  applied  at  r'  and  the  solution  is  sought  at  a 
surface  observation  point  at  range  r.  Then,  it  can  be  shown 
(Ewing,  Jardetsky  and  Press,  1957,  p.40)  that  the  integral 
expression  for  the  vertical  component  of  displacement  at  r 
has  the  form: 

00 

w«,,r,r'>  -  -  jjL-  J  JQ  (kr)J0  <kr')  dk  (2-3) 

"'ll 


and  a, 3  are  the  compressional  and  shear  wave  velocities  re¬ 
spectively.  Using  the  identity 
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(1) 


(2) 


(Z)  = 


*1"  (3, 


(2-5) 


Equation  (2-3)  can  be  expanded  into  a  sum  of  four  integrals 


of  the  form: 
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Now,  by  replacing  k  by  the  complex  argument  k  +  it,  the  path 
of  integration  can  be  deformed  into  the  first  or  fourth 
quadrant  as  shown  in  Figure  4  where  <  denotes  the  Rayleigh 
wave  pole  (i.e./  F(<)  =  0).  Now  for  each  of  the  integrals  in 


Equation  (2-6) ,  the  contour  must  be  closed  in  the  quadrant 
for  which  the  integral  along  the  arc  (i.e.  Tj  or  T  2 )  van¬ 
ishes.  This  choice  is  dictated  by  the  asymptotic  (i.e.  large 
argument)  behavior  of  the  Hankel  functions  which  are  given  by 

»0l)  <*»  *  ViS 

-  VS 

Now,  r  +  r'  >  0  and  consequently  the  first  integral  in  (2-6) 
converges  only  if  the  contour  is  closed  in  the  first  quadrant 
and,  since  there  are  no  singularities  there,  it  makes  no  con¬ 
tribution  to  the  solution.  By  the  same  reasoning,  the  con¬ 
tour  for  the  second  integral  in  (2-6)  is  always  closed  in  the 
fourth  quadrant  yielding  a  residue  contribution  from  the  sim¬ 
ple  Rayleigh  wave  pole  at  k  =  < .  However,  the  evaluation  of 
the  third  and  fourth  integrals  in  Equation  (2-6)  depends  on 
the  sign  of  r-r"  and  thus  is  different  depending  on  whether 
r  >  r'  (i.e.  diverging)  or  r  <  r'  (i.e.  converging).  Thus, 
for  the  diverging  case,  the  third  integral  in  (2-6;  makes  no 
contribution  and  combining  the  residue  contributions  from 
the  pole  at  k  =  <  in  the  second  and  fourth  integrals  gives 

W  (a) ,r , r ')  =  ^-  H^2)  (<r)  Jn(<r')  r>r'  (2-8) 

°  2y  F-(1C)  o  0 

where  F'(<)  denotes  differentiation  with  respect  to  k  fol¬ 
lowed  by  evaluation  at  k.  Similarly,  for  the  converging 
case,  the  fourth  integral  in  (2-6)  does  not  contribute  and 
combining  the  residue  contributions  from  the  second  and 
third  integral  gives 

W  (cu  ,r ,  r ')  =  ~  Jn  (<r)  H^2)(<r')  r  <  r'  (2-9) 

°  F'(<)  U  U 


.  IT 

ix 

e  e 


.  IT 

-iX 

e  e 


(2-7) 


19 


Thus, the  converging  Rayleigh  wave  solution  can  be  obtained  by 
interchanging  the  functional  dependence  of  the  source  and  pro¬ 
pagation  terms  in  the  solution  for  diverging  Rayleigh  waves; 
and  it  follows  that  the  solution  for  converging  Rayleigh 
waves  on  a  multilayered  halfspace  can  be  written  down  directly 
from  Equation  (2-1)  as 


Wj  (w,r) 
Uj  (a>,r) 


-iir  Ar  (oj)  Jq  (k^r)  Sj  (oj ) 
(k.r) 

e.  (<a)  -= — J W. (w,r) 

3  J0(kjr)  3 


(2-10) 


where  r2 

Sj  («)  =  J  p(r',a»)  H^2)  (k.r")  r"  dr".  (2-11) 

rl 

Comparing  Equations  (2-10)  and  (2-1) ,  it  can  be  seen  that  the 
converging  wave  solution  is  characterized  by  a  Bessel  function 
propagation  term  as  opposed  to  the  more  familiar  Hankel  func¬ 
tion  propagation  term  associated  with  the  diverging  wave 
solution.  Thus,  for  the  converging  wave  solution,  the  ampli¬ 
tude  of  the  vertical  component  of  motion  reaches  its  maximum 
at  the  origin  (r=0)  while  the  horizontal  component  vanishes 
at  the  center  point  as  would  be  expected  by  symmetry. 

Clearly,  for  an  observation  point  located  within  the  loaded 
region,  the  total  Rayleigh  wave  displacement  will  be  given  as 
a  superposition  of  contributions  arising  from  the  two  wave 
types  given  by  Equations  (2-1)  and  (2-10) . 


2.3  SCALING  LAWS  FOR  RAYLEIGH  WAVES 

With  reference  to  Equation  (2-1) ,  it  can  be  seen  that 
at  a  fixed  distant  observation  point,  the  characteristics  of 
the  Rayleigh  waves  are  determined  primarily  by  the  behavior 
of  the  source  ( S  ( cu ) )  and  Rayleigh  wave  site  response  (AR(w)) 
functions.  In  this  section,  scaling  laws  will  be  derived  for 
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each  of  these  functions  and  then  combined  to  assess  the  de¬ 
pendence  of  the  Rayleigh  wave  spectrum  on  source  and  site 
parameters . 

2.3.1  Source  Scaling 

In  order  to  define  scaling  laws  for  the  Rayleigh  wave 
source  function,  S(<u),  it  is  first  necessary  to  specify  the 
functional  dependence  of  the  airblast  load.  Following  Brode 
(1968)  this  will  be  approximated  as 


p(r,t)  =  iPg(l-T)  (ae"aT+  be'Sx+  ce*YT; 


(2-12; 


where  a,  b,  c,  a,  S,  y  are  pressure-dependent  constants 
usually  given  in  tabular  form,  APg  is  the  peak  shock  pressure 
given  by 


=  3152W  +  240l£  +  Q>0215 


(2-13) 


for  APg  in  psi,  yield  W  in  megatons  and  range  r  in  kilofeet 
and  t  is  the  reduced  time  defined  by 


t  -  t 


(2-14) 


where  t  is  the  shock  arrival  time  and  is  the  duration  of 
s  p 

the  positive  phase.  It  follows  from  (2-12)  and  (2-14)  that 
the  Fourier  transform  of  the  airblast  load  can  be  written  in 
the  form 


-io)t 


p(r,o)  =  APgDD  e 


s  a 

a+icoD. 


+iojD. 


3+iojD. 
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(2-15) 


Consider  first  the  low-f requency  limit.  Now 
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lim  J.(kr)  — 1  (2-16) 

co-0  u 

and  it  follows  from  Equation  (2-2)  that 


lim  S(co) 
co-0 
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lim  p(r',o>)  r '  dr'. 
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(2-17) 


But,  by  the  definition  of  the  Fourier  transform, 


lim  p(r',cu) 
co-0 

and  thus 


t )  dt 


lim  S(co) 
co-0 


t)  r'  dt  dr' 


(2-18) 


(2-19) 


which  is  easily  recognized  as  the  total  impulse  delivered  by 

the  airblast.  Now,  by  Equation  (2-13)  the  peak  pressure  AP 

is  invariant  under  cube -root  scaling  and,  conseauently ,  at 

}, 

the  same  scaled  distance  the  impulse  is  proportional  to  W3 
(i.e.,  D  )  .  Moreover,  the  area  over  which  it  is  applied  is 
proportional  to  r 2  and  thus  proportional  to  W3.  Thus,  the 
total  impulse  is  directly  proportional  to  yield,  or 


lim  s  (o>)  ~  w  (2-20) 

co-0 

independent  of  the  site  characteristics. 

In  the  high-frequency  limit.  Equation  (2-15)  simpli¬ 
fies  to 


.  .  ,  ,  .  _  -lcut  a+b+c  , «  . 

lim  p(r,co)  =  AP_  e  s  — - -  .  (2-21) 
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But  from  (2-12),  the  condition  p(r,0)  =  APs  implies  a+b+c  = 

1  and  it  follows  that 
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(2-22) 
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Furthermore,  when  to  is  large,  kr  is  also  large  over  the  range 
of  integration  rQ  to  »  and  thus  the  Bessel  function  in  (2-2) 
can  be  replaced  by  its  asymptotic  approximation 


J0(kr)  *  VS  Sin(kr  +  l) 
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where  c  denotes  the  Rayleigh  wave  phase  velocity  (c  =  to/ k)  . 
It  follows  that  in  the  high-frequency  limit  Equation  (2-2) 
takes  the  form 


lim 

Qj  -*  QO 


S(to)  =  - 


r^APg  (r") 


(r" ) 


(r')e 


-10) 


)- 


5  dr " 


(2-24) 


Now  tg,  r "  >  0  and  it  follows  that  the  integrand  in  the  sec¬ 
ond  term  of  (2-24)  is  highly  oscillatory  in  the  high-fre¬ 
quency  limit  and  thus  the  contribution  of  this  term  to  the 
solution  is  negligible.  Likewise,  in  the  first  integral  in 
(2-24),  over  the  range  in  which  the  quantity  t  (r)  -  r7c  is 
rapidly  varying,  this  term  will  contribute  little  to  the 
solution.  It  follows  then  that  in  the  high-frequency  limit, 
the  solution  depends  only  on  the  behavior  of  the  first  inte¬ 
gral  in  the  vicinity  of  r  =  r  where  r  is  the  point  at  which 


fc;(r)  E  vtrr 
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(2-25) 


where  the  prime  denotes  differentiation  with  respect  to  r 
and  V(r)  is  the  propagation  velocity  of  the  airblast  at  r 
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(such  a  point  will  always  exist  in  the  range  of  integration 
Tq '  00  provided  c  >  Vg  where  Vg  is  the  sound  velocity  in  air). 
This  contribution  can  be  evaluated  by  the  principle  of  sta¬ 
tionary  phase  (Papoulis,  1962,  p.140)  to  obtain 


lim 


/• 

j  r'*  APs(r') 

Jro 


-iaj  tg(r') 
e  L 


APs(?) 


-ioi 

e 


t  (r) 
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and  it  follows  from  (2-24)  that 


lim  l  S  ( co ) 


APS(?) 


(2-27) 


Thus,  at  high  frequencies  the  Rayleigh  wave  source  spectrum 

is  proportional  to  <x>~2.  The  yield  and  site  dependence  of  the 

spectral  amplitude  level  at  high-frequencies  can  be  evaluated 

from  (2-27)  once  the  dependence  of  t  on  r  has  been  specified. 

s  k 

For  example,  at  ranges  greater  than  159W3m  (W  in  megatons) 
Brode  (1968)  suggests 

r% 

t  (r)  -  (2-28) 

s 

and  it  follows  that 

IS 

t"  (r)  ~  n-  -  (2-29) 

s 

But  r  ~  which  implies  t" (r )  ~  W  3  and,  using  the  fact 

that  APs  is  invariant  under  cube-root  scaling,  it  follows 
from  (2-27)  that 

lim  [S(w)j~W3.  (2-30) 

0J-*  <30 
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A  similar  line  of  reasoning  can  be  followed  to  deduce 
the  dependence  of  the  high-frequency  spectral  amplitude  level 
on  the  site  geology,  which  enters  through  the  Rayleigh  wave 
phase  velocity  c.  Suppose  that  at  a  fixed  high-frequency  Wq 
and  yield  W  the  phase  velocity  at  one  site  is  c^  and  at  an¬ 
other  site  is  C£  where  $2  =  k  c^.  Mow,  if  the  airblast  pro¬ 
pagation  velocity  V(r)  =  c^  at  r  =  r^,  it  follows  from  (2-28) 
that  V(r)  =  Cj  at  r  =  rj  where  =  r^/k*.  Then  it  follows 
from  (2-27)  that 


and,  since  AP  is  a  decreasing  function  of  r,  lim  |  S  ( w ) '  in- 

S  (j- 

creases  as  c  increases.  Moreover,  in  the  high-frequency 
limit,  c  is  proportional  to  the  near-surface  shear  wave 
velocity  and  thus  the  high-frequency  Rayleigh  wave  spectral 
amplitude  level  will  increase  as  the  near-surface  shear  wave 
velocity  increases. 

These  scaling  laws  are  summarized  graphically  in  Fig¬ 
ures  5  and  6.  Figure  5  shows  a  schematic  representation  of 
the  Rayleigh  wave  source  spectra  corresponding  to  explosive 
yields  W^,  with  W2  >  W^.  By  Equation  (2-19)  the  spectral 
amplitude  levels  approach  constant  values  at  low  frequencies 
(i.e.,  the  respective  total  impulses)  and  this  level  scales  as 
the  first  power  of  the  yield.  At  high-frequencies,  by  Equa¬ 
tion  (2-27),  the  spectra  are  proportional  to  i-:  and  the 

spectral  amplitude  level  at  a  fixed  high-frequency  is  pro- 

v 

portional  to  W 3.  This  implies  the  existence  of  a  corner  fre- 

—  , 

quency  which  scales  like  W  3,  such  as  that  shown  on  tnis 

figure.  Figure  6  shows  a  schematic  representation  of  the 
Rayleigh  wave  source  spectra  corresponding  to  two  different 
sites  characterized  by  high-frequency  phase  velocities  c^,c ^ 


RELATIVE  SOURCE  SPECTRAL  AMPLITUDE 


with  c 2  >  c^.  Since  the  low-frequency  spectral  amplitude 
level  is  independent  of  the  site  geology  and  the  high-fre¬ 
quency  spectral  amplitude  level  increases  with  increasing  c, 
it  follows  that  the  site  with  the  higher  value  of  c  must  be 
characterized  by  a  higher  corner  frequency  leading  to  a  dif¬ 
ference  in  spectral  composition  such  as  that  shown  on  this 
figure. 

2.3.2  Site  Scaling 

The  Rayleigh  wave  site  response  function,  Ar(gj) ,  in¬ 
corporates  all  the  complexity  of  the  subsurface  geology  and, 
consequently,  in  the  general  case  the  relationship  between 
the  spectral  composition  of  AR(co)  and  the  physical  properties 
of  the  site  is  not  a  simple  one.  However,  for  the  special 
case  of  a  single  layer  over  a  half space,  a  simple  scaling  law 
can  be  derived  which  provides  insight  into  the  response  char¬ 
acteristics  of  some  more  complicated  geologic  models,  includ¬ 
ing  those  of  interest  in  MX  siting. 

The  simple-layer-over-halfspace  model  selected  for 
analysis  is  shown  in  Figure  7,  where  it  is  indicated  that  the 
J.ayer  thickness  (H)  and  shear  wave  velocity  in  the  layer  (3^) 
will  be  considered  as  free  parameters  to  be  varied.  Figures 

8  and  9  show  the  fundamental  mode  AR(co)  as  a  function  of  fre¬ 
quency  for  different  values  of  H  and  3^.  Figure  8  shows  the 
variation  with  H  at  a  fixed  value  of  3^  (316  m/sec) ;  Figure 

9  shows  the  variation  with  3^  at  a  fixed  value  of  H  (12  m) . 

In  general,  the  AR  functions  are  characterized  by  a  well  de¬ 
fined  corner  frequency  below  which  the  response  drops  off 
very  rapidly  and  above  which  the  response  increases  at  a  rate 
which  asymptotically  approaches  the  first  power  of  frequency 
(Harkrider  et  al ,  1974).  A  more  important  observation  with 
regard  to  the  present  investigation  is  that  these  corner 
frequencies  correlate  very  well  with  the  frequencies  of  min¬ 
imum  group  velocity  which  are  indicated  by  the  locations  of 
the  tips  of  the  arrows  on  these  two  figures.  The 
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Figure  7.  Geologic  model  for  the  parametric  study  of  the 
Rayleigh  wave  site  response  function. 


FREQUENCY,  Hz 


Figure  8 .  Variation  of  Ar  spectral  response  as  a  function  of 
H  for  a  fixed  value  of 
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Figure  9.  Variation  of  Ar  spectral  response  as  a  function  c 
3.  for  a  fixed  value  of  H. 


significance  of  this  fact  is  that  it  can  be  shown  that  for  a 
wide  range  of  physical  property  choices,  the  frequency  of 
minimum  group  velocity  is  directly  proportional  to  B^  and  in¬ 
versely  proportional  to  H.  In  fact,  for  the  fundamental  mode, 
this  frequency  (f  )  can  be  approximated  very  closely  by  the 
relation  (Mooney  and  Bolt,  1966) 

fc  =  TTTh  '  (2-32) 


While  this  relation  strictly  applies  only  to  the  single -layer- 
over-half  space  problem,  it  has  been  found  that  it  can  be  ap¬ 
plied  to  MX-type  geologic  environments  with  the  following 
modifications.  Let  H  denote  the  depth  to  the  first  signifi¬ 
cant  discontinuity  in  the  shear  wave  velocity  (e.g.,  soil/rock 
interface)  and  6  the  depth-averaged  shear  wave  velocity  given 
by 


where  EH.  =  H.  Then,  the  corner  frequency  of  the  fundamen- 
i  1 

tal  mode  site  response  function  can  be  approximated  as  (Auld 
and  Murphy,  1979) 


f  ~ 
"c  ~ 


(2-34) 


2.3.3  Rayleigh  Wave  Scaling  Relations 

It  now  remains  to  combine  the  scaling  laws  for  S(oj) 
and  A^(cu)  so  that  the  dependence  of  the  total  Rayleigh  wave 
spectrum  on  yield  and  site  characteristics  can  be  assessed. 
Figure  10  shows  a  schematic  representation  of  the  AR  func¬ 
tions  for  two  hypothetical  sites.  It  can  be  seen  that  in 
both  cases  the  amplitude  response  increases  as  oi3  up  to 


32 


Schematic  representation  of  the  Ar  spectral 
response  at  two  hypothetical  sites. 


i 

» 
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widely-separated  corner  frequencies  above  which  the  response 
increases  more  gradually  (i.e.  ^  cu)  .  These  spectra  have  been 
multiplied  by  the  simplified  source  spectra  of  Figure  5  to 
illustrate  the  range  of  possible  Rayleigh  wave  yield  scaling 
conditions.  Figure  11  shows  the  results  of  multiplying  the 
source  functions  of  Figure  5  by  the  '  Rayleigh  wave  site 
response  function  of  Figure  10.  In  this  case,  the  corner  fre¬ 
quency  of  the  site  response  function  lies  below  the  corner 
frequency  of  the  source  functions  and  as  a  result  the  peak 
Rayleigh  wave  spectral  response  scales  with  the  source  func¬ 
tion.  That  is,  the  frequency  of  maximum  spectral  amplitude 

level  scales  as  W  ^  and  the  maximum  amplitude  level  is  di- 

y, 

rectly  proportional  to  W  .  However,  as  is  illustrated  in 
Figure  12,  the  result  is  quite  different  when  these  same  two 
source  functions  are  multiplied  by  the  A^1^  site  response 
function  of  Figure  10.  In  this  case,  the  source  function 
corner  frequencies  lie  below  the  corner  frequency  of  the 
site  response  function  which  now  controls  the  frequency  of 
maximum  Rayleigh  wave  response.  Consequently,  in  this  case, 
the  dominant  Rayleigh  wave  frequency  is  independent  of  yield 
and  the  maximum  spectral  amplitude  is  proportional  to  W^3  . 


Although  the  above  examples  have  contrasted  the  effects  of 
different  site  response  characteristics,  it  should  be  noted 
that  since  the  source  function  corner  frequency  is  propor¬ 
tional  to  W  ^ ,  there  will  be  a  yield  level  for  every  site 
above  which  the  dominant  Rayleigh  wave  frequency  remains 
constant.  It  will  be  shown  in  the  following  section  that 
for  the  range  of  MX  prototype  geologies  considered  to  date, 
this  yield  threshold  appears  to  be  of  the  order  of  a  few 
hundred  kilotons . 


2.4  RAYLEIGH  WAVE  GROUND  MOTION  ESTIMATES  FOR  SOME 
PROTOTYPE  MX  SITES 

Most  of  the  siting  alternatives  currently  being  con¬ 
sidered  for  the  MX  system  involve  deployment  in  one  or  more 
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FREQUENCY  (Hz) 


Figure  11.  Schematic  diagram  illustrating  the 

Rayleigh  wave  yield  scaling  relations 
for  site  )  .  Here  fs  denotes  the 
site  corner  frequency  and  fw^  denotes 
the  source  function  corner  frequency 
for  yield  W , . 
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Figure  12.  Schematic  diagram  illustrating  the 

Rayleigh  wave  yield  scaling  relations 
for  site  Here  fs  denotes  the 

site  corner  frequency  and  fw^  denotes 
the  source  function  corner  frequency 
for  yield  W . . 


» 
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valleys  located  somewhere  in  the  Basin  and  Range  Province  of 
the  southwestern  U.S.  Thus,  it  is  possible  to  specify  the 
range  of  site  geologies  that  may  be  encountered  and  use  these 
site  descriptions  to  estimate  the  characteristics  of  the  low- 
frequency  ground  motion  environment  to  be  expected  under  a 
variety  of  conditions.  Subsurface  geologic  profiles  for  the 
five  sites  selected  for  analysis  by  the  Data  Analysis  Working 
Group  (1978)  are  shown  in  Figures  13  and  14.  It  can  be  seen 
that  they  are  all  characterized  by  a  relatively  thick  surfi- 
cial  layer  of  soil  (90  to  1000  m)  representing  the  alluvium- 
filled  valleys.  The  physical  properties  which  have  been  as¬ 
signed  to  the  various  geologic  units  are  summarized  in  Table 
1.  Note  that  for  these  models  the  first  significant  discon¬ 
tinuity  in  the  shear  wave  velocity  always  occurs  at  the  soil/ 
rock  interface.  Thus,  although  the  water  table  represents  a 
significant  discontinuity  in  compressional  wave  velocity,  it 
corresponds  to  only  a  relatively  minor  change  in  the  shear 
wave  velocity. 

The  fundamental  mode,  Rayleigh  wave  site  response 
functions  for  these  five  sites  are  shown  in  Figures  15 
through  17.  It  can  be  seen  that  the  computed  corner  fre¬ 
quencies  (i.e.,  fc)  vary  over  a  broad  range  (0.3  to  2.3  Hz). 
Moreover,  if  H  is  defined  as  the  total  soil  thickness  and  3 
the  average  soil  velocity,  then  it  can  be  shown  that  the  cor¬ 
ner  frequencies  for  these  five  cases  scale  in  very  good 
agreement  with  Equation  (2-34).  In  terms  of  the  examples 
shown  in  Section  2.3.3  above,  however,  the  question  remains 
as  to  whether  the  dominant  Rayleigh  wave  frequency  will  be 
controlled  by  the  source  or  site  response  characteristics  for 
attack  scenarios  of  interest.  In  order  to  address  this  ques¬ 
tion,  a  series  of  calculations  have  been  performed  assuming 
that  each  of  the  five  sites  was  loaded  with  the  airblast  to 
be  expected  from  a  single,  one  megaton  surface  burst. 
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Figure  17.  Fundamental  mode  Rayleigh  wave  site  response 
function  for  Case  5. 


As  in  the  case  of  the  previous  HE  calculations  (Murphy, 
1978)  the  airblast  load  for  the  Rayleigh  wave  calculations 
has  been  approximated  by  a  series  of  simple  functions  of  the 


form 


-a.r  -6 .  (t-r . ) 

P(r,t)  =  C^e  1  e  x  H(t-ri) 


(2-35) 


in  order  to  simplify  the  numerical  calculations.  Here,  the 
subscript  i  denotes  the  fact  that  a  series  (i.e.  ,  four)  of 
such  functions  is  used  to  describe  the  overpressure  in  various 
distance  ranges,  denotes  the  arrival  time  of  the  airblast 
at  radius  r  and  H  is  the  unit  step  function.  Taking  the 
Fourier  transform  of  Equation  (2-35)  and  substituting  into 
Equation  (2-2)  leads  to  the  following  expression  for  the 
diverging  component  of  the  Rayleigh  wave  source  function: 


Sj  (<u) 


-iojt . 

C^e  1(8i-ioj) 

Sl  S*  +■  co* 


r >e~ar  JQ(kjr")cos  $ 


-r 


r'e  ar  (kjr') sin  ^  r'dr" 


r  "dr ' 


(2-36) 


where  j  denotes  the  mode  number,  is  the  average  airshock 
propagation  velocity  across  the  surface  disk  bounded  by  r^_^ 
and  r^  and  t^  is  the  arrival  time  of  the  shock  at  r^_^.  A 
similar  relation  holds  for  the  converging  component  of  the 
source  function  (i.e.,  Sj  (w)  ,  Equation  (2-11))  with  the  Hankel 
function  replacing  the  Bessel  function  in  the  integrands  of 
Equation  (2-36) .  The  parameters  of  the  overpressure  fit  for 
this  case  have  been  computed  by  scaling  those  from  the  100- 
ton  fit  (Murphy,  1978)  to  a  yield  of  one  megaton  and  are 
given  in  Table  2.  It  can  be  seen  from  this  table  that  the 
innermost  disk  extending  from  the  origin  out  to  246  m  is 
loaded  with  an  overpressure  which  is  of  considerably  lower 
amplitude  than  that  which  would  ordinarily  be  predicted  in 
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this  range  for  an  explosion  of  this  yield.  This  reflects  the 
observation  that  loading  inside  this  scaled  radius  does  not 
seem  to  be  efficiently  coupled  into  radiated  surface  waves. 
The  reason  for  this  effect  is  not  known  at  this  time,  but  may 
be  related  to  energy  losses  associated  with  crater  formation 
and  other  strong  nonlinear  interactions  which  are  known  to 
take  place  in  the  immediate  vicinity  of  ground  zero  (Murphy, 
1978)  . 

The  above  defined  loading  conditions  have  been  used  to 
compute  Rayleigh  wave  ground  motion  time  histories  for  each 
of  the  five  prototype  MX  sites.  The  resulting  fundamental¬ 
mode  displacement  components  at  a  fixed  range  of  1800  m  are 
shown  in  Figures  18  through  26.  For  Cases  1  through  4,  dis¬ 
placement  time  histories  computed  using  finite  difference 
codes  are  also  shown  for  purposes  of  comparison.  The  finite 
difference  solutions  for  Cases  1  through  3,  which  are  shown 
in  Figures  18  through  23,  were  computed  using  a  nonlinear 
finite  difference  code  in  which  che  geologic  models  were 
driven  by  the  complete  airblast  loading  expected  from  a  one 
megaton  surface  burst  (Sandler,  1978) .  The  finite  difference 
solutions  for  Case  4,  on  the  other  hand,  were  computed  using 
a  linear,  elastic  finite  difference  code  in  which  the  geo¬ 
logic  model  was  driven  by  a  modified  airblast  load  comparable 
to  that  defined  by  Equation  (2-35)  and  Table  2  (Auld,  1978) . 
Although  the  details  of  the  modeling  assumptions  are  differ¬ 
ent  enough  to  preclude  absolute  amplitude  comparisons  at  this 
time,  the  general  character  of  the  low-frequency  ground  mo¬ 
tions  ,  including  dominant  frequency  and  relative  displacement 
amplitude  levels  between  sites,  would  be  expected  to  agree 
fairly  well.  This  seems  to  be  true  for  Cases  1  through  3, 
for  which  it  appears  that  the  fundamental  mode  Rayleigh  wave 
dominates  the  displacement  motions  at  this  range.  However, 
the  agreement  between  the  computed  motions  for  Case  4  is  not 
nearly  as  good,  and  this  is  somewhat  puzzling  in  that  the 
airblast  ioading  functions  used  in  the  finite  difference  and 
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alytical  Rayleiyh  wave  and  finite  difference  displace- 
ies,  Case  1,  vertical  component,  W  =  1  Mt,  R  =  1800  m. 


ANALYTICAL  RAYLEIGH  WAVE 


ANALYTICAL  RAYLEIGH  WAVE 


Comparison  of  analytical  Rayleigh  wave  and  finite  difference  displace 
ment  time  histories,  Case  2,  vertical  component,  W  =  1  Mt,  R  =  1800  m 
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analytical  Rayleigh  wave  and  finite  difference  displace 
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igure  22.  Comparison  of  analytical  Rayleigh  wave  and  finite  difference  displace 
ment  time  histories,  Case  3,  vertical  component,  W  =  1  Mt,  R  =  1800  m 
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tories.  Case  4,  vertical  component,  W  =  1  Mt,  R  =  1800  m. 
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arison  of  analytical  Rayleigh  wave  and  finite  difference  displace- 
time  histories,  Case  4,  radial  component,  W  =  1  Mt,  R  =  1800  m. 


analytic  Rayleigh  wave  calculations  were  most  nearly  the  same 
in  this  case.  Thus,  although  the  dominant  frequencies  are 
still  comparable,  the  general  character  of  the  motions  is 
quite  different.  This  may  be  an  indication  that  the  higher 
mode  Rayleigh  waves  are  making  a  significant  contribution  in 
this  case,  although  preliminary  theoretical  calculations  which 
include  the  contributions  from  the  first  higher  mode  do  not 
seem  to  support  this  interpretation. 

In  any  case,  it  can  be  seen  that  the  predominant  fre¬ 
quencies  of  the  Rayleigh  wave  and  finite  difference  solutions 
are  in  good  agreement  with  each  other  and  also  with  the  site 
response  corner  frequencies  shown  in  Figures  15  through  17. 
Thus,  it  can  be  concluded  that  the  corner  frequencies  of  the 
Rayleigh  wave  source  functions  corresponding  to  a  one  megaton 
surface  burst  lie  below  the  site  response  corner  frequencies 
in  all  of  these  cases.  In  fact,  parametric  calculations  in¬ 
dicate  that  this  is  true  for  all  yields  above  a  few  hundred 
kilotons.  That  is,  for  yields  above  this  l^vel  the  predomi¬ 
nant  Rayleigh  wave  frequency  can  be  expected  to  be  indepen¬ 
dent  of  yield  and  the  amplitude  level  is  expected  to  be  pro- 
portional  to  W 3  at  a  fixed  range  for  these  sites.  This  is 
illustrated  in  Figure  27  which  compares  the  fundamental  mode, 
vertical  component  Rayleigh  wave  displacements  at  a  fixed 
range  of  1800  m  computed  for  Case  2  assuming  modified  air- 
blast  loadings  corresponding  to  1  and  30  megaton  surface 
bursts.  It  can  be  seen  that  the  dominant  frequency  of  the 
motion  remains  constant  and  that  the  amplitude  level  in¬ 
creases  as  the  cube  root  of  the  ratio  of  the  yields,  as  pre¬ 
dicted  . 

In  summary,  available  evidence  suggests  that  the 
dominant  frequency  and  relative  amplitude  level  of  the  low- 
frequency  ground  motions  to  be  expected  at  a  variety  of  pro¬ 
totype  MX  sites  from  megaton  level  surface  bursts  can  be 
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Comparison  of  vertical  component  Rayleigh  wave 
displacement  time  histories  for  yields  of  1  and 
30  MT ,  Case  2,  R  =  1800  m. 


reliably  estimated  using  elastic  Rayleigh  wave  theory.  In 
particular,  the  dominant  frequency  can  be  approximated  as 
S/2H,  where  H  denotes  the  total  soil  thickness  at  the  site 
and  3  is  the  average  shear  wave  velocity  of  the  soil  column. 
For  yields  above  a  few  hundred  kilotons ,  this  dominant  fre¬ 
quency  will  be  independent  of  yield  and  the  amplitude  of  the 
Rayleigh  wave  ground  motion  at  a  fixed  range  will  be  propor¬ 
tional  to  W1^  . 

The  results  of  the  above  analysis  suggest  a  possible 
simplified  means  of  ranking  candidate  MX  sites  according  to 
the  expected  severity  of  the  low-frequency  ground  motion  en¬ 
vironment.  Thus,  for  each  site  under  consideration,  S(fc) 
can  be  approximated  from  (2-27)  for  the  yield  of  interest 
(for  W  greater  than  about  200  kt)  and  multiplied  times  A„(f  ) 
to  obtain  an  estimate  of  the  relative  displacement  amplitude 
level  to  be  expected  at  the  predominant  Rayleigh  wave  fre¬ 
quency  of  the  site  (i.e.  fc)  .  Moreover,  since  the  computed 
time  histories  are  approximately  sinusoidal,  the  relative 
displacement  amplitude  levels  could  then  be  converted  to 
estimates  of  the  corresponding  relative  particle  velocity 
amplitude  levels  by  simply  multiplying  by  the  quantity  2Tf  . 

2.5  MULTIBURST  SURFACE  WAVE  EFFECTS  ON  MISERS  BLUFF  II-2 

The  Misers  Bluff  II  (MBII)  experiments  were  a  series 
of  two  high  explosive  surface  burst  tests  conducted  by  the 
Defense  Nuclear  Agency  in  a  small  scale  valley  structure  at 
Planet  Ranch  in  Arizona.  The  second  test  conducted  at  this 
site  (MBII-2)  consisted  of  an  array  of  six  simultaneous  120 
ton  ANFO  explosions  deployed  in  a  hexagonal  pattern.  Figure 
28  shows  a  schematic  drawing  of  the  source  array  together 
with  some  representative  receiver  locations  as  well  as  a  plane 
layered  approximation  to  the  subsurface  geology  at  the  site. 

It  can  be  seen  that  the  site  is  characterized  by  a  70  m  thick 
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57  m 
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p  =  2.4  gm/cm"* 

igure  28.  Source/receiver  pattern  and  sub¬ 
surface  geologic  profile  for  the 
Misers  Bluff  II-2  experiment. 


layer  of  alluvium  (saturated  below  the  water  table  at  a 
depth  of  13  m)  overlying  basalt.  An  unusual  feature  of  the 
ground  motion  recordings  taken  interior  to  this  source  array 
is  the  presence  of  a  large  burst  of  energy  arriving  at  the 
center  of  the  array  about  800  milliseconds  after  detonation 
time.  This  arrival  is  illustrated  in  Figure  29  which  shows 
the  observed  vertical  displacement  records  measured  interior 
to  the  source  array  along  a  common  azimuth.  It  can  be  seen 
that  the  maximum  peak-to-peak  displacement  at  the  center  of 
the  source  array  (i.e.  r  =  0)  is  about  five  times  larger 
than  that  observed  at  the  station  located  midway  between  the 
origin  and  the  source  radius  (i.e.  r  =  50  m) .  Analysis  of 
the  complete  set  of  data  recorded  interior  to  the  source 
array  indicates  that  the  diagnostic  characteristics  of  this 
signal  are:  (i)  that  it  is  propagating  inward  toward  the 
center  of  the  array  with  a  phase  velocity  of  around  500  m/sec, 
(ii)  that  the  amplitude  of  the  vertical  component  of  the 
motion  increases  to  a  maximum  at  the  center  of  the  array, 
as  opposed  to  the  amplitude  of  the  horizontal  component  of 
motion  which  reaches  a  minimum  at  that  point  and  (iii)  that 
the  signal  is  approximately  axisymmetric .  Since  no  such 
arrival  was  observed  on  the  companion  single  burst  experi¬ 
ment  conducted  at  this  site  (MBII-1)  and  since  linear  super¬ 
position  of  the  motions  from  the  six  explosions  is  unable 
to  explain  the  observed  spatial  variation  of  the  displacement 
peaks,  it  has  been  concluded  that  this  arrival  was  generated 
by  the  nonlinear  interaction  between  the  multiple  sources. 
Moreover,  the  observed  phase  velocity  and  character  of  the 
signal  suggest  that  it  may  be  a  surface  wave  and  consequently 
the  arrival  has  been  modeled  as  a  Rayleigh  wave  converging 
on  the  origin.  Now,  despite  the  fact  that  the  airblast 
loading  is  known  to  be  very  complicated  in  this  case,  the 
symmetric  nature  of  the  late-time  arrival  of  interest 
here  suggests  that  its  source  may  be  axisymmetric  to 
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a  first  order  of  approximation.  Consequently,  the  MBII-2 
observations  have  been  theoretically  simulated  by  loading  a 
surface  ring  of  the  geologic  model  shown  in  Figure  28  with  a 
simple,  impulsive  (i.e.,  a  =  10  sec  ^)  overpressure  approxi¬ 
mation  of  the  form 


p(r,t)  =  pQ  e“a(t_t0)  H(t-t0) 


(2-37) 


where  tg  is  a  free  parameter  defining  the  time  required  for 
the  interacting  shock  waves  from  the  six  events  to  form  an 
approximately  axisymmetric  loading  condition.  In  this  case, 
the  source  integral  of  Equation  (2-11)  has  the  form 


S  (oj)  = 


a+icu 


kr')  r'  dr' 


(2-38) 


which  can  be  integrated  analytically  to  obtain 


S(ai) 


„  -iwtrt 
P0e  0 

(a+iu)  5c 


ip0e-iait0 
(a+icu)  5c 


£r2  J1(kr2)  -  r1  J1(kr1)j 

J^r2  Y1(kr2)  -  r1  Y1(kr1)j.  (2-39) 


The  converging  fundamental  mode  Rayleigh  waves  produced  by 
applying  three  different  ring  loads  centered  on  the  source 
radius  (100  m)  to  the  geologic  model  of  Figure  28  are  shown 
in  Figures  30  through  32.  It  can  be  seen  that  the  theoreti¬ 
cal  prediction  for  the  case  shown  in  Figure  32  (i.e.,  r^  =  90 
m,  r2  =  110  m,  tQ  *  0.4  sec)  fits  the  observations  very  well. 
In  particular,  the  phase  velocity,  dominant  frequency  and 
surface  wave  waveforms  are  nearly  exactly  duplicated  while 
the  predicted  relative  amplitude  levels  are  within  the  scat¬ 
ter  in  the  experimental  data. 


62 


VERTICAL  DISPLACEMENT 


VERTICAL  DISPLACEMENT,  cm 


OBSERVED 


THEORETICAL  RAYLEIGH  WAVE 


R  =  12.5  m 
(AV  2209) 


R  =  25  m 
(AV  2138) 


4.0 


0 


R  =  50  m 
(AV  2153) 


-4.0 


l— I  1  1  1  !  1  1  1  1  1 

0  0.4  0.8  1.2  1.6  2.0 

TIME ,  sec 


Figure  31.  Comparison  of  observed  MBII-2 

vertical  displacement  data  with 
theoretical  converging  Rayleigh 
wave  solution;  distance  measured 
from  center  of  source  array  ( t q  = 

0.4  seconds,  r,  =  75  m,  r7  =  125  m) . 
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Figure  32.  Comparison  of  observed  MBII-2  vertical 
displacement  data  with  theoretical  con¬ 
verging  Rayleigh  wave  solution;  distance 
measured  from  center  of  source  array 
(tn  =0.4  seconds,  r,  =  90  m,  r_  =  110  m) . 


These  results  provide  strong  evidence  that  the  300 
millisecond  arrival  observed  on  MBII-2  is  a  converging  Rayleigh 
wave  produced  by  an  axisymmetric  surface  load  acting  over  a 
relatively  small  area  centered  on  the  source  locations.  The 
question  still  remains,  however,  of  how  this  axisymmetric 
loading  is  set  up  by  the  interacting  shock  waves  from  the  six 
MBII-2  sources.  For  example,  in  order  to  match  the  observed 
arrival  times  in  Figure  32,  it  was  necessary  to  assume  a 
value  of  tg  =  0.4  seconds.  In  view  of  the  fact  that  the  air- 
blast  reaches  the  center  of  the  source  array  in  less  than  0.1 
seconds,  this  is  a  very  late  time  phenomena.  It  is  inter¬ 
esting  to  note  that  a  major  spall  closure  did  occur  at  about 
0.4  seconds  on  the  companion  single  explosive  event  conducted 
at  this  site  (MBII-1)  suggesting  that  the  interacting  spall 
closures  from  the  multiple  explosions  may  be  the  source  of 
this  arrival.  However,  more  detailed  analyses  of  the  mea¬ 
sured  MBII-2  near-source  data  will  be  required  in  order  to 
provide  a  definitive  test  of  this  hypothesis. 

In  summary,  multiburst  loadings  resulting  from  at 
least  some  charge  spacings  can  produce  long-period  surface 
wave  motions  which  seem  to  have  no  counterparts  on  single 
burst  recordings.  This  suggests  that  the  possible  influence 
of  such  phenomena  on  the  ground  motion  environment  for  MX 
should  be  carefully  evaluated  on  a  case-by-case  basis. 
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Ill .  VALLEY  REVERBERATION 


3 . 1  INTRODUCTION 

It  has  been  recognized  from  the  early  planning  stages 
that  the  geologic  environment  of  potential  MX  sites  (i.e.,  al¬ 
luvium-filled  valleys)  is  quite  different  from  that  encoun¬ 
tered  in  the  Minuteman  siting.  As  a  result,  questions  have 
been  raised  concerning  the  issue  of  whether  the  ground  motion 
environment  in  the  valley  produced  by  near-surface  explosions 
will  be  significantly  altered,  and  in  particular  enhanced,  as 
a  result  of  the  presence  of  the  valley  boundaries.  For  ex¬ 
ample,  Figure  33  shows  a  schematic  geologic  cross  section  of  a 
typical  Basin  and  Range  valley.  It  can  be  seen  that  at  the 
edges  of  the  valley  the  alluvium  is  terminated  by  boundary 
faults,  giving  rise  to  a  significant  seismic  impedance  con¬ 
trast  at  the  steeply-dipping  alluvium/rock  interface.  Thus, 
it  can  be  expected  that  at  least  some  of  the  outgoing  seismic 
energy  induced  by  a  surface  explosion  will  be  reflected  back 
from  the  valley  boundaries  (particularly  for  the  horizontally 
propagating  surface  waves)  and  contribute  to  the  ground 
motion  experienced  at  a  target  point  in  the  valley.  This  is 
in  contrast  to  the  situation  at  the  Minuteman  sites  where  the 
subsurface  geologic  structure  is  more  nearly  laterally  homo¬ 
geneous  and,  consequently,  contributions  to  the  ground  motion 
due  to  reflections  from  vertical  boundaries  would  be  expected 
to  be  negligible. 

An  example  of  the  effects  of  valley  boundaries  is  pro¬ 
vided  by  Figure  34  where  the  radial  component  seismograms  re¬ 
corded  at  a  station  in  the  Las  Vegas  Valley  (SE-6)  and  at  a 
station  outside  the  valley  (PAH)  at  about  the  same  distance 
from  the  same  Nevada  Test  Site  (NTS)  explosion  are  compared. 
It  can  be  seen  that  the  upper  trace  which  was  recorded 
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Schematic  geologic  cross  section  for  a  Basin  and  Range  valley  structure. 


Comparison  of  seismocjrams  recorded  at  about  the  same  epicentral  distance 
inside  (SE-6)  and  outside  (PAH)  the  Las  Vegas  valley. 


within  the  valley  is  much  more  complex  and  of  significantly 
longer  duration  than  the  lower  trace  which  was  recorded  out¬ 
side  the  valley.  The  inference  is  that  multipathing  (i.e., 
multiple  reflections)  of  energy  along  the  valley  boundaries 
is  the  source  of  this  dramatic  increase  in  complexity  and 
duration. 

During  the  past  several  years  a  number  of  specific  ex¬ 
perimental  and  theoretical  investigations  have  been  conducted 
with  the  objective  of  defining  the  effect  of  valley  bound¬ 
aries  on  the  ground  motion  environment  to  be  expected  at  tar¬ 
get  points  inside  the  valley.  The  ultimate  goal  of  these 
studies  is  to  develop  a  better  understanding  of  the  phenome¬ 
nology  so  that  an  approximate  synthesis  model  can  be  formu¬ 
lated  which  can  be  used  in  predicting  the  environment  in  the 
MX  valleys  of  interest. 

This  chapter,  first,  summarizes  the  results  of  prior 
investigations  which  suggest  that  alterations  in  the  ground 
motion  can  be  anticipated  in  such  environments  and,  then, 
focuses  on  analyses  of  ground  motion  measurements  from  the 
Misers  Bluff  II  (MB-II)  explosions  and  their  implications 
with  regard  to  the  effects  of  valley  boundaries. 

3.2  RESULTS  FROM  PREVIOUS  STUDIES 

Prior  investigations  of  ground  motion  near  the  bound¬ 
aries  of  alluvial  valleys  have  included  both  observational 
and  theoretical  studies;  we  will  describe  here  the  pertinent 
results  from  both  types  of  study.  Experimental  studies  to 
date  have  centered  on  the  analysis  of  ground  motion  data  re¬ 
corded  in  valleys  from  NTS  underground  explosions  as  well  as 
data  recorded  on  specially  designed  arrays  of  seismometers 
deployed  to  monitor  the  Misers  Bluff  II  (MB-II)  explosions. 

The  MB-II  studies  will  be  discussed  more  completely  in  Section 

3.3  below.  Although  valuable  information  has  been  obtained 
from  both  these  data  sources,  it  should  be  noted  that 
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neither  are  directly  applicable  to  the  definition  of  the  MX 
environment.  That  is,  the  NTS  data  were  measured  from  con¬ 
tained  underground  explosions  while  the  MB^II  data  were  mea¬ 
sured  from  small-scale  surface  HE  detonations  conducted  in  a 
"miniature  MX  valley".  Therefore,  the  primary  goal  of  our 
consideration  of  these  data  has  been  to  develop  a  better 
understanding  of  the  phenomenology  so  that  extrapolations  can 
be  made  to  the  full-scale  environment  of  interest. 

The  first  data  bearing  on  this  problem  was  alluded  to 
earlier  and  consisted  of  ground  motion  recordings  made  in  Las 
Vegas  from  contained  underground  explosions  located  on  NTS 
some  150  km  away  and  outside  of  the  Las  Vegas  Valley  struc¬ 
ture.  The  Las  Vegas  Valley  is  a  prominent  northwest- trending 
topographic  depression  situated  in  the  southern  part  of  the 
Basin  and  Range  Province.  The  city  of  Las  Vegas  is  located 
in  the  southern  portion  of  this  elongated  valley  in  the  cen¬ 
ter  of  a  rectangular  basin  which  is  about  30  km  wide  and  40 
km  long.  The  city  has  been  well  instrumented  to  record 
ground  motions  from  underground  nuclear  tests  conducted  at 
NTS.  Figure  35  shows  the  locations  of  a  selected  subset  of 
these  stations  and  Figure  36  shows  the  radial  component  veloc¬ 
ity  time  histories  recorded  at  these  stations  from  the 
CARPETBAG  explosion.  These  records  are  typical  and  show  pro¬ 
nounced  evidence  of  interference  at  late  arrival  times  con¬ 
sistent  with  what  would  be  expected  from  the  simultaneous 
arrival  of  energy  from  several  different  directions  (Murphy 
and  Hewlett,  1975) . 

In  order  to  determine  the  direction  of  approach  and 
relative  amplitude  of  these  secondary  arrivals,  the  ground 
motion  time  histories  of  Figure  36  were  velocity  filtered 
using  the  beam-steer  process  which  is  routinely  used  in  the 
detection  and  analysis  of  signals  recorded  at  large  seismic 
arrays  such  as  the  Montana  Large  Aperture  Seismic  Array  (LASA) 
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The  concept  underlying  this  process  is  quite  simple.  Suppose 
an  arrival  of  interest  (e.g.,  a  P  wave)  is  traveling  across 
the  array  from  some  unknown  direction  with  phase  velocity  c. 
Now,  for  each  possible  direction  of  approach,  8  ,  there  will 
be  a  unique  set  of  delays,  At. (6)  describing  the  propagation 
across  the  station  locations  at  velocity  c.  Thus,  if  the 
waveforms  recorded  at  each  station  are  time-shifted  by  an 
amount  At^(0)  and  added  together,  then  energy  arriving  from 
direction  0  with  velocity  c  will  add  coherently  while,  in 
general,  energy  traveling  in  other  directions  or  at  other 
velocities  will  add  incoherently  and  tend  to  cancel.  Thus, 
the  value  of  0  for  which  a  given  arrival  achieves  maximum 
amplitude  is  associated  with  the  approach  azimuth  of  the  com¬ 
ponent  of  that  arrival  traveling  with  velocity  c.  The  re¬ 
sponse  of  such  arrays  is  discussed  more  fully  in  Section  3.3 
below. 


In  the  Las  Vegas  Valley  study,  since  visual  analyses 
of  the  seismograms  suggested  that  the  reflected  arrivals  were 
predominantly  Rayleigh  waves  with  periods  ranging  from  only 
about  2.0  to  4.0  seconds,  the  simple  beam-steer  process  de¬ 
scribed  above  was  modified  to  take  advantage  of  this  addition¬ 
al  information  (Murphy  and  O'Brien,  1978).  Thus,  since  the 
arrivals  of  interest  were  confined  to  a  relatively  narrow 
frequency  band,  the  measured  data  were  filtered  through  an 
eighth-order  bandpass  filter  with  a  transfer  function  of  the 
form  , 


T  (S)  = 


u)  S 

c 


S*+  a>cS 


+  0)  ‘ 


3-1) 


where  S  is  the  Laplace  transform  variable  and  =  2-Tf^  with 
f  =  0.3  Hz.  This  effectively  eliminated  shorter  period  body 
wave  phases  which  were  not  of  primary  interest  in  this  case. 
Next,  in  order  to  isolate  the  contribution  to  the  observed 
ground  motion  due  to  Rayleigh  waves  arriving  from  a  selected 
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approach  azimuth,  8,  use  was  made  of  the  fact  that  the  Ray¬ 
leigh  wave  particle  motion  is  restricted  to  the  radial-verti¬ 
cal  plane  oriented  along  the  approach  azimuth.  That  is,  the 
radial  component  of  motion  with  respect  to  9  was  computed  by 
applying  a  rotation  of  axes  to  the  measured  orthogonal,  hori¬ 
zontal  components  of  motion.  Finally,  the  filtered  and  ro¬ 
tated  data  from  the  array  stations  were  summed  with  inter¬ 
station  time  delays  appropriate  for  a  plane  wave  approaching 
from  direction  8  at  the  Rayleigh  wave  velocity  corresponding 
to  a  frequency  of  0.3  Hz. 

The  results  of  applying  this  filtering  procedure  to 
the  CARPETBAG  data  of  Figure  36  are  shown  in  Figure  37.  In 
this  figure,  zero  degrees  indicates  the  azimuth  back  toward 
the  explosion  source  and  the  angle  is  measured  clockwise  from 
north.  It  can  be  seen  that  the  output  of  the  filter  in  this 
case  clearly  indicated  some  strong  secondary  bursts  of  energy 
arriving  from  directions  other  than  along  the  azimuth  back  to 
the  source.  The  strongest  of  these  was  an  arrival  origina¬ 
ting  at  an  azimuth  between  160  and  180  degrees  which  propa¬ 
gated  back  up  the  line  toward  NTS.  The  general  amplitude 
level  associated  with  this  arrival  was  found  to  be  comparable 
to  that  of  the  direct  Rayleigh  waves,  suggesting  that  the 
reflection  coefficient  at  the  boundary  which  is  the  source  of 
this  phase  must  be  close  to  unity  for  incident  Rayleigh  waves 
with  periods  near  3.0  seconds.  The  only  other  prominent 
secondary  burst  of  energy  which  could  be  easily  identified 
(c.f..  Figure  37)  was  a  still  later  arrival  originating  at  an 
azimuth  of  about  120  degrees  which  was  quite  similar  in 
character  to  that  originating  at  180  degrees. 

In  order  to  complete  the  identification,  it  was  nec¬ 
essary  to  correlate  these  arrivals  with  the  structural  bound¬ 
aries  of  the  valley.  Figure  38  is  a  map  of  the  Las  Vegas 
Valley  showing  the  locations  of  the  surrounding  mountain 
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ranges  with  respect  to  the  instrumentation  array  in  Las  Vegas. 
It  can  be  seen  that  the  valley  is  terminated  on  the  southeast 
by  the  McCullough  mountain  range,  and  the  relative  arrival 
time  of  the  secondary  signal  from  an  azimuth  of  180  degrees 
was  found  to  be  consistent  with  it  being  a  reflection  from 
this  boundary  if  the  average  Rayleigh  wave  propagation  veloc¬ 
ity  at  a  period  of  3.0  seconds  is  assumed.  Perhaps  even  more 
interesting  is  the  observation  that  the  arrival  time  of  the 
phase  propagating  across  the  array  from  an  azimuth  of  120 
degrees  was  consistent  with  it  being  a  multiply  reflected 
arrival  from  Frenchman  Mountain  following  the  rather  compli¬ 
cated  path  shown  on  this  figure.  If  this  latter  identifica¬ 
tion  is  correct  then  this  arrival  must  correspond  to  a  fairly 
localized  focusing  of  energy,  since  its  amplitude  seems  to  be 
too  large  to  be  consistent  with  the  hypothesized  non-planar 
reflection. 

The  Las  Vegas  data  analysis  presented  above  clearly 
illustrates  the  fact  that  the  presence  of  the  valley  bound¬ 
aries  can  significantly  modify  the  ground  motion  experienced 
in  the  valley.  However,  it  is  not  directly  relevant  to  the 
MX  valley  reverberation  question  in  that  the  source  was  lo¬ 
cated  some  150  km  away  outside  of  the  valley.  In  order  to 
examine  the  problem  on  a  more  appropriate  scale,  data  re¬ 
corded  from  the  NTS  event  CALABASH  in  Yucca  Valley  have  also 
been  analyzed.  Figure  39  shows  the  locations  of  CALABASH  and 
the  stations  that  were  deployed  to  measure  the  near-field 
ground  motions  from  this  event.  It  can  be  seen  from  this 
figure  that  the  instrumentation  array  encompassed  stations 
located  at  about  the  same  distance  (7  to  8  km)  both  outside 
and  inside  the  valley.  Figure  40  shows  a  vertical  subsurface 
geologic  section  along  the  line  A-A'  on  Figure  39  which  indi¬ 
cates  that  the  alluvium  in  the  area  gradually  thins  to  the 
west  until  it  terminates  against  the  underlying  tuff  forma¬ 
tion  at  the  valley  boundary.  Figure  41  shows  a  comparison  of 
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Vertical  geologic  section  for  the  CALABASH  event. 
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the  radial  component  velocity  time  histories  measured  for 
CALABASH  at  stations  752  (in  the  valley)  and  755  (on  the 
valley  boundary) .  It  can  be  seen  that  although  these  stations 
were  located  at  about  the  same  distance  from  the  source,  there 
was  a  significant  difference  in  the  effective  duration  with 
the  closer,  in-valley  station  (752)  showing  clear  evidence  of 
later  secondary  arrivals.  This  difference  is  even  more  ap¬ 
parent  when  the  envelopes  of  the  signals  in  a  narrow  fre¬ 
quency  band  are  compared.  Figure  42  shows  a  block  diagram  of 
the  envelope  detector  which  was  selected  to  process  the  re¬ 
corded  data.  The  passband  around  1  Hz  was  chosen  for  analy¬ 
sis  in  this  case  on  the  basis  of  a  visual  identification  of 
the  characteristic  frequencies  of  the  late-time  motions  re¬ 
corded  on  the  array.  Figure  43  shows  the  1  Hz  envelopes  of 
the  signals  shown  in  Figure  41,  and  Figure  44  shows  the  en¬ 
velopes  for  another  pair  of  stations  one  of  which  was  inside 
(746)  and  one  of  which  was  outside  (756)  the  valley.  In  both 
cases  there  is  clear  evidence  of  strong  secondary  arrivals  at 
the  stations  inside  the  valley  which  appear  to  have  no  coun¬ 
terparts  on  the  records  of  adjacent  stations  located  outside 
the  valley.  This  strongly  suggests  that  a  significant  amount 
of  energy  is  being  reflected  back  from  the  valley  boundary 
and  modifying  the  motion  inside  the  valley. 

We  next  turn  to  the  theoretical  studies;  we  will  re¬ 
turn  to  observational  results  from  the  MB-II  studies  in 
Section  3.3.  The  theoretical  investigation  of  valley  rever¬ 
beration  is  a  difficult  problem  because  of  the  complexity  of 
the  geologic  structure  involved.  For  example,  it  is  question¬ 
able  whether  valleys  such  as  those  discussed  above  can  be 
adequately  modeled  without  going  to  a  full  three  dimensional 
description  of  the  valley  boundaries.  Although  such  a  sim¬ 
ulation  is  probably  within  the  current  state  of  the  art,  it 
would  undoubtedly  be  very  costly  and  time  consuming,  parti¬ 
cularly  in  view  of  the  fact  that  there  is  so  little  previous 


Block  diagram  of  the  bandpass  filter  envelope  circuit. 


Figure  44.  Comparison  of  the  1  Hz  envelope  for  stations  inside  (746)  and  outside 
(756)  Yucca  Valley. 


experience  at  modeling  structures  of  this  complexity  at  a 
realistic  scale.  For  these  reasons,  the  two  studies  conducted 
to  date  (Orphal  and  Hancock,  1978;  Cherry,  et  al. ,  1978)  have 
consisted  of  finite  difference  simulations  of  the  response  of 
highly  simplified,  two-dimensional  approximations  to  a  valley 
structure.  The  results  of  these  preliminary  studies  are  sum¬ 
marized  in  the  following  paragraphs. 

Orphal  and  Hancock  (1978)  considered  the  response  of  a 
20  km  wide  elliptical  alluvial  valley  of  infinite  length 
assuming  plane-strain  symmetry.  The  geometry  and  loading 
conditions  for  the  problem  are  illustrated  in  Figure  45  which 
shows  that  the  entire  valley  was  simultaneously  loaded  with 
an  overpressure  pulse  of  constant  peak  amplitude  (i.e.,  1  bar) 
and  a  pulse  shape  consistent  with  that  expected  at  that  peak 
pressure  from  a  1  megaton  surface  burst  (Brode,  1968) .  This 
loading  was  taken  as  a  crude  approximation  to  the  effect  of 
a  multiburst  attack  on  an  MX  valley.  The  simulation  was 
carried  out  using  the  Physics  International  PISCES  2D  ELK 
continuum  mechanics,  finite  difference  computer  code.  The 
problem  was  further  simplified  by  assuming  that  the  response 
of  the  geologic  materials  to  the  loading  was  perfectly  elas¬ 
tic  everywhere.  Figure  46  shows  the  computed  horizontal 
particle  velocity  time  histories  as  a  function  of  range  from 
the  center  of  the  valley  (i.e.,  r  =  0) .  Since  the  companion 
simulation  run  in  the  absence  of  the  valley  boundaries  has 
not  yet  been  performed,  it  is  difficult  to  make  a  quantita¬ 
tive  assessment  of  how  much  the  valley  boundaries  are  affec¬ 
ting  the  ground  motion  response.  However,  several  strong, 
low-frequency  secondary  phases  can  be  identified  which  are 
clearly  propagating  back  in  toward  the  center  of  the  valley. 
This  is  more  clearly  illustrated  in  Figure  47  which  shows 
"snapshots"  of  the  horizontal  particle  velocity  as  a  function 
of  position  in  the  valley  at  progressively  later  times.  It 
can  be  seen  that  as  time  advances  a  standing-wave  pattern 


Figure  45.  Geometry  and  loading  conditions  for  elliptical  valley  approximation. 


Figure  46.  Computed  horizontal  particle  velocity  seismograms 
as  a  function  of  distance  from  the  valley  center. 
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seems  to  be  building  up  in  the  valley,  although  the  problem 
was  not  run  long  enough  to  allow  the  pattern  to  fully  develop 
In  any  case,  it  seems  clear  that  a  significant  amount  of 
energy  is  being  reflected  back  from  the  valley  boundaries, 
leading  to  an  increase  in  the  duration  of  motion  which  is  at 
least  qualitatively  consistent  with  the  observed  data. 

Cherry,  et  al.  (1978)  considered  a  cylindrical  valley 
approximation  and  used  a  special  purpose,  elastic  finite 
difference  code  to  compute  the  ground  motion  response  both 
with  and  without  a  valley  boundary.  In  both  cases,  the  sur¬ 
face  area  outside  the  expected  final  crater  region  (i.e., 
region  of  strong  nonlinear  response)  was  loaded  with  a  Brode 
(1968)  airblast  appropriate  for  a  single,  megaton  yield  sur¬ 
face  burst.  Figure  48  shows  the  subsurface  geological  model 
used  in  these  simulations.  It  can  be  seen  that  this  model 
was  more  complicated  than  that  used  by  Orphal  and  Hancock 
(c.f.,  Figure  45)  and  was,  in  fact,  chosen  to  be  a  plane 
layered  representation  of  a  typical  Basin  and  Range  valley 
structure.  An  initial  simulation  was  made  in  which  these 
layers  were  taken  to  have  infinite  horizontal  extent  (i.e., 
no  valley  boundaries) ,  and  this  served  as  a  basis  for  com¬ 
parison.  In  the  second  simulation,  a  vertical  boundary  was 
inserted  at  a  range  of  3  km  from  the  center  of  the  valley, 
and  it  was  surrounded  by  the  underlying  basement  rock  (i.e., 
granite) .  Figures  49  through  54  show  comparisons  of  the 
vertical  and  horizontal  components  of  the  particle  velocity 
at  three  different  ranges  (i.e.,  0.9,  2.03  and  2.96  km  from 
the  center  of  the  valley)  computed  with  and  without  the 
valley  boundary.  Again,  it  can  be  seen  that  the  valley  bound¬ 
aries  significantly  affected  the  ground  motion  response  in 
the  valley,  particularly  with  regard  to  the  duration  of  the 
longer-period  surface  waves.  For  example,  the  upper  traces 
in  Figures  49  and  50  indicate  that  at  a  range  of  0.9  km  the 
effective  duration  of  motion  was  less  than  10  seconds  in  the 
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absence  of  the  valley  boundary.  However,  at  that  same  loca¬ 
tion  for  the  valley  simulation  (i.e.,  the  lower  traces  in 
these  figures)  the  amplitude  of  the  lower  frequency  (~0.6 
Hz)  surface  wave  motion  20  seconds  after  detonation  time, 
when  the  calculation  stopped,  was  still  comparable  to  that  of 
the  direct-arriving  surface  wave  energy.  This  indicates  that 
a  substantial  amount  of  energy  is  being  reflected  back  from 
this  vertical  discontinuity  in  physical  properties.  In  gen¬ 
eral,  the  effect  of  the  valley  boundaries  on  the  peak  veloc¬ 
ity  values  were  seen  to  be  much  less  pronounced  than  the 
effect  on  duration.  However,  an  exception  to  this  occurred 
for  the  vertical  velocity  at  the  station  very  close  to  the 
boundary  (i.e.,  r  ■  2.96  km).  Here,  as  can  be  seen  from 
Figure  54,  the  peak  velocity  for  the  problem  with  the  valley 
boundary  was  actually  about  a  factor  of  two  smaller  than  that 
which  was  observed  in  the  absence  of  the  valley  boundary. 

This  effect  was  found  to  be  due  to  destructive  interference 
between  the  incident  and  reflected  Rayleigh  waves  and  has 
been  used  by  Cherry,  et  al.  (1978)  to  infer  a  Rayleigh  wave 
reflection  coefficient  of  about  0.5  for  this  geologic  model. 

3.3  MISERS  BLUFF  II  REVERBERATION  STUDIES 
3.3.1  Data  Analyses 

As  noted  above,  the  Misers  Bluff  II  (MBII)  explosions 
•'ere  a  series  of  two  HE  tests  designed  to  simulate  the  blast 
ground  motion  environment  in  an  MX  valley.  One  element  of 
the  MBII  study  was  specifically  designed  to  test  the  effects 
on  ground  motion  of  reverberation  of  seismic  waves  in  the 
"miniature  MX  valley"  at  Planet  Ranch  in  Arizona.  Instrumen¬ 
tation  for  this  element  of  the  study  consisted  of  the  seisric 
stations  shown  in  Figure  55.  The  stations  for  the  MBII-1 
explosion  are  shown  as  sclid  circles  and  the  stations  for  the 
MB I I- 2  explosion  are  shown  as  open  circles;  the  station  at 
the  southern  edge  of  the  valley  was  in  place  for  both 
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Figure  55.  Bedrock  contour  map  for  Misers  Bluff  II  site  with 
locations  for  individual  explosions  (+)  and  re¬ 
cording  stations  (•)  for  the  MBII-1  single  burst 
event  and  (o)  for  the  MBII-2  multiple  burst  event. 
Contour  lines  give  bedrock  elevations  in  feet; 
local  average  surface  elevation  is  670  feet. 
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explosions.  Each  station  included  two  horizontal  instruments 
designed  to  respone  to  radial  (i.e.,  along  the  shot/station 
line)  and  transverse  (i.e.,  perpendicular  to  the  shot/station 
line)  motions.  In  addition,  four  stations  for  each  explo¬ 
sion  were  equipped  with  vertical  instruments.  The  station 
near  the  southern  edge  of  the  valley  also  included  three- 
component  down-hole  instrumentation.  Velocity  records  from 
the  broad-band  instruments  were  integrated  to  obtain  dis¬ 
placement  time  histories  which  were  used  in  our  analysis. 

Also  shown  in  Figure  55  are  bedrock  contours  of  the 
valley.  The  contours  on  this  figure  show  the  bedrock  eleva¬ 
tion  in  feet;  the  average  surface  elevation  at  the  site  is 
about  670  feet.  It  can  be  seen  that  even  in  the  center  of 
the  valley  the  alluvium  thickness  varies  between  20  and  70 
meters.  In  particular,  note  that  there  is  a  local  high  in 
the  bedrock  topography  located  between  the  source  and  the 
station  array.  As  might  be  expected,  this  local  anomaly, 
which  was  mapped  after  the  design  of  the  reverberation  array 
had  been  finalized,  further  complicates  the  interpretation 
of  an  already  complex  data  set. 

Figure  56  shows  a  comparison  of  ground  motion  records 
measured  at  a  range  of  about  300  m  from  two  different  100  ton 
HE  experiments.  The  lower  trace  was  measured  from  the  Pre- 
Mine  Throw  IV  (PMT  IV)  test  conducted  at  the  Nevada  Test  Site 
and  has  an  effective  duration  of  only  about  one  second.  By 
contrast,  the  upper  trace  recorded  at  about  the  same  range 
from  the  MBII-1  test  in  the  Planet  Ranch  valley  shows  an  ef¬ 
fective  duration  of  at  least  ten  seconds.  Thus,  there  is 
clear  evidence  of  reverberation  of  seismic  energy  within  the 
valley  structure.  Another  observation  that  can  be  made  from 
this  recording  is  that  the  characteristic  frequency  associ¬ 
ated  with  the  secondary  arrivals  is  significantly  lower  than 
that  of  the  initial,  direct  arrivals.  This  is  illustrated 
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more  clearly  in  Figure  57  which  shows  the  result  of  passing 
one  of  the  recorded  MBII-1  displacement  time  histories 
through  eighth-order  bandpass  filters  centered  at  4.0  and  2.0 
Hz.  It  can  be  seen  that  most  of  the  4.0  Hz  energy  is  concen¬ 
trated  near  the  beginning  of  the  signal,  while  the  2.0  Hz 
component  of  the  motion  maintains  a  relatively  continuous 
level  over  the  entire  10  seconds  of  recording.  It  seems 
likely  that  this  difference  in  characteristic  frequency  re¬ 
flects  differences  in  average  properties  of  the  associated 
propagation  paths.  This  is  illustrated  in  Figure  58  which 
shows  the  fundamental  Rayleigh  wave  site  response  functions 
for  the  Planet  Ranch  site.  Now,  the  direct  arrivals  travel 
along  a  path  for  which  the  average  thickness  of  alluvium  is 
around  40  m,  giving  a  characteristic  Rayleigh  wave  frequency 
of  4.0  Hz.  The  late-time  reverberation  signals,  on  the  other 
hand,  may  traverse  portions  of  the  valley  where  the  alluvial 
thickness  is  70  m  or  more,  leading  to  a  characteristic 
Rayleigh  wave  frequency  closer  to  2.0  Hz. 

We  processed  the  MBII  array  data  using  a  procedure  of 
frequency  filtering  and  beam-steering  like  that  used  in  the 
analysis  of  the  Las  Vegas  array  data  described  above.  A 
schematic  describing  the  beam-steering  process,  as  used  here, 
is  shown  in  Figure  59.  Stations  in  the  recording  array,  rep¬ 
resented  by  the  black  circles,  are  uniformly  separated  by  a 
distance,  d.  The  wave  motion  recorded  on  the  array  is  as¬ 
sumed  to  consist  of  plane  waves  which  travel  at  a  velocity, 
c,  and  propagate  at  an  azimuth  angle  9  relative  to  the  line 
of  stations.  Assuming  that  there  is  no  change  in  waveform 
across  the  array  and  the  motion  at  station  1  is  represented 
by  F  ( t) ,  then  the  motions  at  other  stations  in  the  array  are 
simply  delayed  in  time  by  an  amount  equal  to  the  time  for  the 
wavefront  to  travel  from  the  position  intersecting  station  1 
to  a  position  intersecting  the  other  stations.  Then,  the 
wave  motion  at  stations  in  the  array  can  be  represented  as: 


f  *2.0  Hz 


Time,  seconds 


Figure  57.  Comparison  of  eighth  order  bandpass  filtered  dis 
olacements  at  4.0  and  2.0  Hz;  MBII-1,  r  =  305  m. 
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inear,  uniformly-spaced  array 


f  x  (t)  =  F  ( t)  ,  f2(t)  =  F  (t-At)  ,  .  . 


(3-2) 


fn(t)  =  F  (t-  [n-1]  At) 

where  At  *  d  cos9/c. 

In  the  beam-forming  process,  recorded  time  histories 
from  individual  stations  in  the  array  are  summed;  however, 
prior  to  summing  time  delays  or  advances  may  be  introduced 
to  "steer"  the  array.  Thus,  for  the  array  shown  in  Figure 
59 ,  advancing  the  time  history  recorded  at  station  n  by  an 
amount  (n-l)d  cos9/c  and  then  summing  causes  the  signals  to 
be  added  in  phase  and  thus  "steers"  the  beam  in  the  direction 
6  from  which  the  signal  arrives.  In  general,  a  time  advance 
t  =  (n-l)d  cos0q/c  for  the  linear  array  can  be  thought  of  as 
steering  the  beam  in  the  direction  Sq  for  signals  of  velocity, 
c . 

A  better  understanding  of  how  the  process  works  can  be 
obtained  by  considering  the  response  of  the  array  to  harmonic 
wave  motion  arriving  from  direction  9.  We  will  assume  a  rime 
advance  given  by  t  =  (n-1) At,  where  At  ®  d  cos0Q/c,  is  in¬ 
troduced  at  station  n  before  the  beam  sum  is  formed.  Taking 
the  amplitude  of  the  harmonic  wave  as  A,  the  beam  sum  is: 

t  ,  Aai“t  .  AeiN<e*AT>-«ltf+_ 

+  Aei  {w  (t+  [N-1]  At)  -<u  [N-1]  At} 

=  Aeiwt  e1  n(a-3)  (3-3) 

n=0 

where  N  is  the  total  number  of  stations  in  the  array  over 
which  the  summation  is  performed,  a  =  jAt  =  ^d  cos  5q/c 
and  3  h  uAt  =  cod  coss/c.  This  sum  can  be  rewritten  so  that: 

,  -  a  <<*?(», i-SMij  \  .  „.4) 

I  sm  ((a-3)/2!  ( 
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The  term  in  braces,  Ssin'T(a'-sl /  '  corresPon<^s  to  t^ie 

ratio  of  the  amplitude  of  the  beam  sum  to  the  signal  ampli¬ 
tude  at  a  single  station  while  the  term  [(a-3 )  / 2]  cor_ 

responds  to  the  change  in  phase  of  the  beam  sum  relative  to 
the  signal  at  the  first  detector. 

In  the  case  of  the  MBII  array,  some  additional  com¬ 
plexity  is  introduced  into  the  problem  by  the  fact  that  we 
are  dealing  with  horizontal  components  of  motion.  The  prob¬ 
lem  is  simplified  by  considering  signals  for  which  the  wave 
particle  motion  is  in  the  plane  of  propagation  separate  from 
signals  for  which  the  particle  motion  is  perpendicular  to  the 
plane  of  propagation.  For  particle  motion  in  the  plane  of 
propagation  (which  is  of  prime  interest  in  the  current  study), 
the  along-the-beam  response  for  an  array  steered  in  the  di¬ 
rection  9q  can  be  written  as: 

^  „  /Q  ,  sin  [N  (ct-S )  /  2]  „i(N-l)  f(a-B ) /2]  „iot  ,,_ex 

*L  -  W  sin  [ia4i/;S]  e  e  (3  51 

where 

AL(90>  =  AR  cos60  "  AT  sin0O  (3-6) 

where  AR  and  are  the  radial  and  transverse  components  of 
the  recorded  ground  motion.  (The  term  Al(0q)  Present  -*-n 
the  amplitude  expression  at  each  station  and  has  been  fac¬ 
tored  out  of  the  sum  to  obtain  Equation  (3-5).)  The  radial 
component  of  motion,  A_,  is  along  the  line  of  stations  and 
the  transverse  component,  A^,  is  perpendicular  to  that  line, 
so  that 


Ar  =  Ar  cos6 
At  *  -AH  sine 


(3-7) 


where  A„  is  the  amplitude  of  the  horizontal  component  in  the 
n 

plane  of  propagation  for  wave  motion  propagating  in  direction 
e  -  i.e.,  Ar  is  the  resultant  of  AR  and  A^ .  Substituting  in 
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Equation  (3-6) , 


w  = 


Ah  COS0  cos9  q  +  Ah  sin0 
Ah  cos  (3-6g) 


sin0 , 


(3-8) 


So : 


L 


*  _,a  fl  x  (sin  [n  (a-8 )  /  2 1  |  _i(N-l)  f(a-3)/2l  iot 

Ah  cos  0-3  q)  ,-sin  !  e  e  * 

(3-9) 


Then  the  ratio  of  the  amplitude  of  the  beam  sum  to  the  ampli¬ 
tude  of  the  horizontal  component  at  a  single  station  for  mo¬ 
tion  in  the  plane  of  propagation  is  just  the  product  of  the 
original  directivity  function  for  the  array  (the  term  in 
braces)  and  a  factor  which  projects  the  motion  into  the  di¬ 
rection  in  which  the  beam  is  steered.  Similar  expressions 
can  be  derived  for  the  perpendicular  response  of  the  beam 
sum  to  signals  with  particle  motion  in  the  plane  of  propaga¬ 
tion  and  for  the  responses  of  the  beam  sum  to  signals  having 
transverse  particle  motion. 


Figure  60  shows  the  predicted  response  of  the  MBII-1 
array  to  harmonic  motion  with  a  frequency  of  4  Hz  and  a 
velocity  of  600  m/sec  for  particle  motion  polarized  in  the 
direction  of  propagation  of  a  signal  from  ?  =  O’.  The  re¬ 
sponse  is  shown  as  a  function  of  the  direction  in  which  the 
beam  is  steered,  9g.  The  response  is  seen  to  consist  of  a 
bread  main  lobe  and  some  lesser  side  lobes.  In  this  case, 
the  response  is  symmetric  about  the  line  of  stations;  for 
signals  propagating  in  other  directions  the  symmetry  normally 
associated  with  the  directivity  of  a  linear  array  is  lost  be¬ 
cause  of  the  factor  projecting  the  motion  onto  the  steering 
direction,  as  will  be  seen  below. 


In  applying  the  beam-steering  process  to  the  MBII-1 
array,  the  recorded  signals  were  not  harmonic  but  consisted 
of  transients  of,  at  most,  a  few  cycles.  Applying  an  eighth- 
order  narrow  band-pass  filter  to  the  motion  makes  the  time 
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histories  appear  more  nearly  harmonic  and,  thus,  should  make 
the  observed  response  of  the  MBII-1  array  more  closely  match 
the  theoretical  response  shown  in  Figure  60.  Similar  proce¬ 
dures  were  used  in  studying  the  reverberations  of  seismic 
waves  in  Las  Vegas  Valley  as  noted  above. 

For  MBII-1  filter  center  frequencies  of  4  Hz  and  2  Hz 
were  selected.  The  phase  velocities  used  to  compute  delay 
times  between  stations  were  varied  over'  a  range  from  about 
0.2  km/sec  to  1.4  km/sec  which  includes  the  range  of  phase 
velocities  consistent  with  dispersion  curves  calculated  for 
the  theoretical  model  of  the  MBII  area  shown  above  in  Figure 
26.  Figures  61  and  62  show  the  results  of  applying  this 
process  to  the  MBII-1  data  for  a  frequency  of  4  Hz  and  a 
phase  velocity  of  600  m/sec  (the  velocity  which  maximized 
the  beam  amplitude) .  The  time  histories  shown  represent  mo¬ 
tions  along  the  direction  in  which  the  beam  is  steered;  1.5 
seconds  of  "dead  time"  has  been  introduced  at  the  start  and 
finish  of  each  time  history  to  permit  the  required  shifting 
of  time  histories  in  the  beam-forming  process.  It  can  be 
seen  that  the  combined  effect  of  the  velocity  and  frequency 
filters  effectively  suppresses  the  late-time  reverberations 
and  isolates  the  direct  arriving  4  Hz  component  of  the  motion. 
It  is  interesting  to  note,  however,  that  the  maximum  array 
response  occurs  at  an  azimuth  (measured  clockwise  from  the 
source  azimuth)  of  between  20  and  40  degrees.  This  can  also 
be  seen  in  Figure  60  where  we  compare  the  observed  response 
pattern  of  the  MBII-1  array  with  its  theoretical  response  to 
harmonic  motion  arriving  from  9  =  O’.  The  observed  data 
points  correspond  to  amplitudes  of  the  peak  motion  on  each 
record  which  have  been  normalized  by  dividing  bv  the  maxi¬ 
mum  array  response.  Clearly,  the  observed  peak  motion 
does  not  appear  to  be  arriving  along  the  shot/station  line. 
Figure  63  shows  a  similar  comparison  for  harmonic  motion  ar¬ 
riving  from  3  =  301 .  The  observed  pattern  in  peak  amplitudes 
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Response  of  the  MB I I - 1  array  to  4  Hz  signals  with 
phase  velocity  c  =  600  m/sec.  Solid  line  shows 
theoretical  response  to  harmonic  waves  arriving 
from  6  *  0°.  Circles  indicate  observed  response 
to  the  filter,  beam-steering  process  with  obser¬ 
vations  normalized  to  the  maximum  predicted  re¬ 
sponse  . 
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Figure  61.  MBII-1  filtered  beam  c  =  600  m/sec,  f  =  4.0  Hz 
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Response  of  the  MBII-1  array  to  4  Hz  signals  with 
phase  velocity  c  =  600  m/sec.  Solid  line  shows 
theoretical  response  to  harmonic  waves  arriving 
from  0  =  30°.  Circles  indicate  observed  response 
to  the  filter,  beam-steering  process  with  obser¬ 
vations  normalized  to  the  maximum  predicted  re¬ 
sponse  . 
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of  the  steered  MBII-1  array  data  is  found  to  match  the  theo¬ 
retical  pattern  much  more  closely  in  the  latter  case. 

Two  possible  explanations  for  this  puzzling  observa¬ 
tion  have  been  considered.  The  first  is  that  the  primary 
surface  wave  energy  has  been  ducted  around  the  previously 
referenced  high  point  in  the  bedrock  topography  (c.f.,  Figure 
55)  and  is,  in  fact,  approaching  from  well  off  the  source 
azimuth.  The  other  possibility  is  that  the  effect  is  an 
artifact  of  the  data  processing  resulting  from  the  fact  that 
the  observed  transverse  component  of  motion  is  not  zero,  as 
assumed  in  the  formulation  of  the  simple  model.  That  is,  it 
has  been  shown  (Murphy  and  Bennett,  1979)  that  adding  an  in- 
phase  transverse  component  of  motion  to  the  direct  arriving 
Rayleigh  wave  signal  can  lead  to  a  directional  response  very 
similar  to  that  observed  in  Figures  61  and  62.  It  is  not 
clear  at  the  present  time  which  of  these  interpretations  is 
correct.  We  have  conducted  some  preliminary  particle  motion 
studies  in  an  attempt  to  better  discern  the  nature  of  this 
motion.  Figure  64  shows  particle  motion  diagrams  for  one 
second  of  time  of  the  signal  causing  the  peak  motion  as  re¬ 
corded  at  station  280  filtered  through  a  4  Hz  band-pass  fil¬ 
ter.  The  motions  are  shown  for  a  horizontal  (R-T)  plane  and 
for  a  vertical  (R-Z)  plane.  The  orientation  of  the  R-Z  plane 
was  varied  between  an  azimuth  of  -20’  and  +605  (relative  to 
the  shot/station  line)  in  20°  steps  to  produce  the  diagrams. 
We  see  from  the  diagrams  in  the  R-T  plane  that  the  transverse 
component  of  the  overall  motion  is  minimized  (though  it  is 
still  fairly  large)  when  the  R-Z  plane  is  at  an  azimuth  be¬ 
tween  20’  and  40°  .  The  particle  motion  in  the  R-Z  plane  at 
this  azimuth  is  seen  to  consist  of  prograde  and  retrograde 
elliptical  motion  probably  associated  with  Rayleigh  waves. 

The  corresponding  particle  motion  in  the  R-Z  plane  for  an 
azimuth  of  O'  appears  more  complex  and  is  less  easily  identi¬ 
fiable  as  Rayleigh  type  motion.  This  tends  to  favor  the  firs 
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explanation  for  the  observed  behavior  (i.e.  ,  the  "direct"  sig¬ 
nal  is  arriving  from  well  off  the  source  azimuth) .  We  per¬ 
formed  similar  analyses  on  several  later  time  segments  from 
individual  records  and  also  beam-summed  time  histories  and 
found  elliptical  particle  motions  typical  of  Rayleigh  waves. 

Given  that  the  direct  arrival  has  been  isolated,  it 
remains  to  demonstrate  that  it  is  a  Rayleigh  wave.  For  this 
purpose,  the  MBII-1  event  has  been  simulated  by  computing  the 
Rayleigh  wave  component  of  the  ground  motion  induced  by  the 
airblast  acting  on  the  surface  outside  the  region  of  strong 
nonlinear  response  (i.e.  ,  the  crater)  (Murphy,  1979).  For  this 
purpose,  the  subsurface  geologic  model  shown  in  Figure  58  was 
used  with  H  =  40  m. 

Figure  65  shows  a  comparison  of  the  MBII-1  observed 
filtered  beam  output  with  the  theoretical,  radial  component, 
fundamental  mode  Rayleigh  wave  displacement  at  a  range  of 
300  m.  The  agreement  between  the  computed  and  observed  mo¬ 
tion  is  excellent,  confirming  that  the  observed  motion  is 
predominantly  an  airblast-induced ,  fundamencal  mode  Rayleigh 
wave . 

Due  to  the  fact  that  this  array  was  designed  primarily 
for  4  Hz,  600  m/sec  phase  velocity  signals,  it  is  not  well 
suited  for  separating  out  the  late-arriving,  lower  frequency 
signals.  For  this  reason,  the  higher  velocity  (1000  m/sec), 
lower  frequency  (2  Hz)  beam  was  less  successful  in  decompo¬ 
sing  the  recorded  motion.  This  is  illustrated  in  Figure  66 
which  shows  the  filter  outputs  over  the  azimuth  range  from  0 
to  180  degrees.  Several  secondary  bursts  of  energy  are  evi¬ 
dent,  but  they  are  spread  out  over  such  a  wide  range  in  azi¬ 
muth  that  a  definitive  identic ication  of  their  origin  is  net 
possible.  This  lack  of  resolution  is  primarily  due  to  the 
fact  that  the  array  aperture  is  too  small  with  respect  to  the 
characteristic  wavelengths  of  the  secondary  signal,  as  shown 
by  Murphy  and  3ennett  (1979)  using  synthetic  seismograms. 
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theoretical  and  observed  displacements;  MBI1-1,  H  -  300  m 
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Considering  this  cause  for  the  lack  of  resolution  of 
the  lower-frequency,  higher-velocity  seismic  waves  found  in 
the  MBII-1  data,  somewhat  better  results  could  be  anticipated 
for  the  MBII-2  explosion  because  of  the  larger  array  aperture 
(150  m  versus  100  m  for  MBII-1).  To  test  the  frequency  con¬ 
tent  of  the  signals  from  the  MBII-2  explosion,  we  processed 
the  ground  motion  recorded  at  a  range  of  420  m  using  a  series 
of  eighth-order,  narrow  band-pass  filters.  The  center  fre¬ 
quencies  were  varied  between  1  Hz  and  5  Hz.  The  results  are 
shown  in  Figures  67  through  69.  Although  some  lower -fre¬ 
quency  energy  can  be  seen  on  the  records,  the  signal  is  gen¬ 
erally  concentrated  in  a  fairly  narrow  frequency  band  sur¬ 
rounding  3  Hz.  Furthermore,  any  signals  occurring  after 
about  the  first  4  to  5  seconds  of  record  appear  to  be  signi¬ 
ficantly  lower  in  amplitude.  The  only  exception  seems  to  be 
a  2-3  Hz  signal  which  shows  up  mainly  on  the  tangential  re¬ 
cord  about  7  seconds  after  shot  time;  and  even  this  is  a 
factor  of  6  lower  in  amplitude  than  the  direct  motion.  This 
clustering  of  signals  near  the  start  of  the  record  and  in  a 
rather  narrow  range  of  frequencies  is  somewhat  different  from 
similar  computations  for  MBII-1  (c.f..  Figure  57)  and  makes 
signal  resolution  particularly  difficult. 

We  applied  the  same  beam-steering,  frequency-filtering 
procedure  to  the  MBII-2  array  data  as  that  used  for  MBII-1. 
Center  frequencies  were  again  set  at  2  Hz  and  4  Hz  and  the 
velocities  used  to  compute  delay  times  were  varied  over  a 
range  appropriate  to  the  velocity  structure  of  the  MBII  site. 
Figures  70  and  71  show  the  results  of  applying  the  beam¬ 
steering  process  for  a  filter  frequency  of  4  Hz  and  a  phase 
velocity  of  600  m/sec;  the  records  show  the  along-the-beam 
component  of  motion.  In  these  figures  one  second  of  "dead 
time"  has  been  added  at  the  start  and  finish  of  each  time 
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Figure  67.  Comparison  of  eighth-order  bandpass  filtered 

radial  displacements  at  center  frequencies,  fc, 
of  1.0,  2.0,  2.5,  3.0,  4.0  and  5.0  Hz;  MBII-2, 

R  =  420  m. 
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history  to  allow  for  the  shifting  involved  in  the  process. 

As  expected,  prominent  signals  cluster  in  the  first  half  of 
the  record;  nevertheless,  a  few  distinct  phases  seem  to  be 
identifiable.  The  dominant  signal  along  the  shot/station 
azimuth,  which  corresponds  to  the  direct  Rayleigh  wave,  is 
again  maximized  at  a  velocity  near  600  m/sec.  In  addition, 
a  rather  prominent  phase  (delayed  about  1  second  relative  to 
the  direct  Rayleigh  wave)  appears  to  be  traversing  the  array 
from  an  azimuth  of  180°.  This  is  also  apparent  in  the  re¬ 
sponse  pattern  of  the  array  shown  in  Figure  72.  Comparing 
the  observed  peak  response  of  the  MBII-2  array  (approximated 
in  this  case  as  a  six  element,  equally-spaced  array)  for  4  Hz 
waves  arriving  from  6  =  0° ,  we  see  a  close  match  out  to  3q  = 
70°  beyond  which  waves  traveling  back  up  the  array  cause  di¬ 
vergence  from  the  predicted  pattern.  The  travel  time  of  this 
phase,  about  1.4  seconds  relative  to  shot  detonation  (after 
removing  a  delay  introduced  by  the  filter) ,  coupled  with  a 
propagation  velocity  near  600  m/sec  indicate  a  total  travel 
path  length  of  about  800-900  m.  This  path  length  closely 
matches  that  for  the  direct  Rayleigh  wave  reflected  at  the 
southern  boundary  of  the  valley  back  to  the  array. 

Figures  73  and  74  show  the  results  of  beam-steering 
and  filtering  the  MBII-2  data  for  a  filter  center  frequency 
of  2  Hz.  The  figures  show  the  along-the-beam  array  response 
for  azimuths  from  0°  to  180°  and  phase  velocities  of  600 
m/sec  and  1000  m/sec.  In  this  case,  we  note  the  presence  of 
a  relatively  prolonged  phase  shortly  after  the  direct  Rayleigh 
wave  which  appears  to  be  maximized  near  an  azimuth  of  60°. 

This  phase  was  not  apparent  for  the  4  Hz  filtered  data  shown 
previously  (c.f.,  Figure  70)  and,  considering  its  travel  time 
and  direction  of  approach,  appears  to  be  associated  with  Ray¬ 
leigh  waves  reflected  at  the  northeastern  and  eastern  bound¬ 
aries  of  the  valley. 
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Figure  72. 


Response  of  the  MBII-2  array  to  4  Hz  signals  with 
phase  velocity  c  =  600  m/sec.  Solid  line  shows 
theoretical  response  to  harmonic  waves  arriving 
from  -3  =  0°.  Circles  indicate  observed  response 
to  the  filter,  beam-steered  process  with  observa¬ 
tions  normalized  to  the  maximum  oredicted  response 


For  the  MBII-2  data  we  again  had  difficulty  resolving 
the  direction  of  approach  for  many  of  the  phases  apparent  on 
the  records.  The  phases  are  generally  seen  to  be  spread  over 
a  fairly  large  range  in  azimuth.  The  larger  aperture  of  the 
MBII-2  array  appears  to  only  marginally  increase  this  resolu¬ 
tion.  This  can  be  seen  by  comparing  the  theoretical  harmonic 
responses  for  the  MBII-1  and  MBII-2  arrays  in  Figures  60  and 
72 •  While  the  peak  response  for  the  MBII-2  array  is  somewhat 
higher  and  its  side  lobes  are  somewhat  smaller,  nulls  occur 
in  nearly  the  same  positions  and  the  main  lobe  is  still  fairly 
broad . 

As  a  final  element  in  processing  the  MBII-2  data,  we 
considered  the  signals  recorded  at  the  station  near  the 
southern  edge  of  the  valley.  Figures  75  and  76  show  the  re¬ 
sults  of  passing  the  records  from  the  uphole  instruments 
through  a  series  of  band-pass  filters  centered  at  frequencies 
between  1  and  5  Hz.  By  comparing  these  figures  to  similar 
figures  for  the  array  station  at  a  range  of  420  m  (c.f.,  Fig¬ 
ures  67  through  69)  several  interesting  features  can  be  noted. 
The  most  obvious  difference  between  the  two  sets  of  records 
is  the  change  in  frequency  content;  the  records  from  the 
station  near  the  valley  edge  show  little  sign  of  the  2-3  Hz 
energy  so  prominent  in  the  array  records.  This  may  be  asso¬ 
ciated  with  changes  in  the  Rayleigh  wave  response  of  the 
instrument  sites  or  possibly  interference  between  direct  and 
reflected  Rayleigh  waves  occurring  at  the  valley-edge  station. 
One  other  feature  of  the  records  in  Figures  75  and  76  is  the 
presence  of  significant  ground  motion  prior  to  the  arrival  of 
the  direct  Rayleigh  wave.  The  mechanism  for  the  generation 
of  this  wave  motion  is  unclear  but  appears  to  have  a  source 
at  the  valley  boundary.  Though  the  amplitude  of  this  motion 
is  somewhat  less  at  all  frequencies  than  that  of  the  direct 
Rayleigh  wave,  it  may  still  be  significant  because  of  the  in¬ 
creased  duration  of  ground  motion  which  it  apparently  produces. 
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Figure  75.  Comparison  of  eighth-order  bandpass  filtered 

radial  displacements  at  center  frequencies,  fc, 
of  1.0,  2.0,  2.5,  3.0,  4.0  and  5.0  Hz;  MBII-2, 

R  =  650  m. 
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Figure  76.  Comparison  of  eighth-order  bandpass  filtered 

vertical  displacements  at  center  frequencies,  fc, 
of  1.0,  2.0,  2.5,  3.0,  4.0  and  5.0  Hz;  MBII-2, 

R  =  650  m. 


131 


3.3.2  Discussion  of  Results 

Of  the  results  presented  here  for  the  Misers  Bluff  II 
valley  reverberation  experiment,  one  which  appears  most  signi 
ficant  to  MX  siting  considerations  is  the  prolonged  duration 
of  ground  motion  observed  at  some  stations  for  limited  fre¬ 
quency  ranges.  The  most  obvious  example  of  this  is  the  2  Hz 
filtered  gound  motion  time  history  for  the  MBII-1  explosion 
shown  in  Figure  57.  The  amplitude  of  the  2  Hz  motion  on  this 
record  at  a  time  10  seconds  after  detonation  is  still  almost 
half  as  large  as  the  maximum  motion  on  the  record.  Such  pro¬ 
longed  duration  is  also  apparent  on  the  beam-steered  data 
from  the  MBII-1  experiment  shown  in  Figure  66  for  a  2  Hz  fil¬ 
ter  frequency  and  phase  velocity  of  1000  m/sec.  In  this  case 
we  again  see  that  motion  10  seconds  after  the  shot  detonation 
are  more  than  one-third  as  large  as  the  peak  motion. 

Although  some  difficulty  was  encountered  in  trying  to 
resolve  late-arriving  signals  because  of  the  relatively  broad 
array  response  and  the  numerous  phases  which  appear  to  be 
arriving,  some  commentary  on  the  nature  of  these  late  ar¬ 
rivals  and  their  propagation  paths  does  seem  warranted.  Firs 
particle  motion  diagrams  prepared  for  the  time  windows  asso¬ 
ciated  with  several  of  the  prominent,  late-time  phases  in 
Figure  66  indicate  elliptical  particle  motion  in  the  R-Z 
plane.  Thus,  we  can  conclude  that  many  of  these  late  ar¬ 
rivals  are  probably  associated  with  Rayleigh  waves.  Secondly 
the  relatively  late  time  of  these  phases  (~  10  seconds)  and 
their  associated  propagation  velocities  (near  1000  m/sec) 
indicate  a  propagation  path  length  which  could  be  as  great  as 
10,000  m.  This  is  much  greater  than  the  dimensions  of  the 
valley  itself  and  indicates  either  reflection  from  boundaries 
outside  the  valley  structure  or,  more  likely,  multiple  re¬ 
flections  within  the  valley  structure.  Assuming  multiple 
reflections  are  responsible  for  the  motion,  reflections  a: 
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as  many  as  10  boundaries  (or  more)  may  be  needed  to  produce 
the  long  path  length.  This  would  imply  an  extremely  efficient 
propagation  path  with  highly  reflective  boundaries.  Further¬ 
more,  if  multiple  reflections  within  the  valley  structure  are 
producing  the  late  time  (~10  seconds)  arrival,  many  of  the 
earlier  phases  adding  to  the  complexity  of  the  beamed  records 
in  Figure  66  are  probably  associated  with  shorter  segments  of 
the  long  multiple  reflection  path.  The  nearly  constant  ampli¬ 
tude  level  associated  with  these  phases  (after  the  first  few 
seconds  of  record)  is  also  indicative  of  propagation  with 
little  loss  of  seismic  energy. 

Another  interesting  feature  of  the  ground  motion  dura¬ 
tion,  which  is  not  completely  understood  at  this  time,  is 
that  the  amplitude  level  at  late  time  relative  to  the  ampli¬ 
tude  of  the  direct  Rayleigh  wave  motion  is  greater  for  the 
single  burst  MBII-1  event  than  for  the  multiple  explosion, 

MB I I- 2 .  This  can  be  seen  most  clearly  by  comparing  the  2  Hz 
filtered,  beam-steered  data  for  the  two  events  (i .  e.,  compare 
Figure  66  with  Figures  73  and  74).  For  the  MBII-1  explosion 
the  late  time  ( — 10—12  seconds  in  Figure  56)  amplitude  level 
is  about  one-third  as  large  as  the  amplitude  level  of  the 
direct  Rayleigh  wave  motion.  In  contrast,  for  the  MBII-2 
explosions  the  amplitude  levels  at  times  after  about  the 
first  5-6  seconds  of  record  in  Figures  73  and  74  are  only 
about  one-eighth  the  amplitude  of  the  direct  Rayleigh  wave 
motion.  In  fact,  the  filtered,  beam-steered  traces  from  the 
MBII-2  explosions  generally  appear  less  complex  than  those 
for  the  MBII-1  explosion.  This  is  really  quite  surprising 
considering  that  the  multiple  explosion  would  intuitively 
seem  to  be  a  more  complex  source.  For  whatever  reason,  the 
multiple  source  (MBII-2)  has  apparently  not  excited  the 
highly  efficient  wave  guide  to  the  same  degree  as  the  single 
source  (MBII-1) .  The  most  obvious  cause  for  such  differences 
would  appear  to  be  the  difference  in  the  locations  of  the 
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events  with  respect  to  the  complex  geologic  structure  of  the 
valley . 

Finally,  one  other  observation  described  above  could 
have  significant  bearing  on  the  problem  of  ground  motion 
duration  in  an  MX  valley.  For  the  station  near  the  valley 
edge,  we  noted  excitation  of  fairly  large  amplitudes  prior  to 
the  arrival  of  the  direct  Rayleigh  wave  (c.f.,  Figures  75  and 
76) .  The  source  of  these  motions  appears  to  be  at  the  valley 
boundary  and  may  result  from  conversion  of  body-wave  to  sur¬ 
face-wave  energy,  or  similar  phenomena,  occurring  there.  The 
implications  of  this  observation  with  regard  to  sites  near 
valley  boundaries  appears  to  be  that  the  generation  of  such 
motion  may  prolong  the  duration  of  significant  levels  of  sur¬ 
face  motion  through  the  body-wave  time  window  and  into  the 
time  where  the  direct  surface  waves  and  subsequent  arrivals 
begin  to  dominate  the  motion.  The  conversion  of  seismic  wave 
energy  at  such  boundaries,  which  appears  to  be  the  source  of 
the  motion,  may  also  have  implications  with  regard  to  the 
reciprocal  problem  of  ground  motion  in  the  valley  due  to 
sources  acting  in  close  proximity  to  its  boundaries.  Such 
conversions  might  be  expected  to  enhance  the  surface  wave 
ground  motions. 
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IV.  SUMMARY  AND  CONCLUSIONS 


4 . 1  SUMMARY 

The  investigations  summarized  in  this  report  have  cen¬ 
tered  on  an  attempt  to  develop  an  improved,  quantitative 
understanding  of  the  late- time,  low-frequency  ground  motion 
environment  for  the  MX  system.  This  effort  has  encompassed 
studies  of  site  and  yield  scaling  laws  applicable  to  this 
component  of  the  ground  motion  as  well  as  an  assessment  of 
the  possibility  that  the  ground  motion  environment  may  be 
significantly  modified  as  a  result  of  siting  the  system  in 
one  or  a  series  of  valleys. 

A  detailed,  theoretical  analysis  of  the  characteris¬ 
tics  of  Rayleigh  waves  produced  by  near-surface  explosions 
was  presented  in  Section  II,  where  scaling  laws  were  derived 
which  describe  the  dependence  of  the  airblast-induced,  Rayleigh 
wave  motion  on  both  the  yield  of  the  explosion  and  the  site 
geology.  These  concepts  were  then  applied  to  an  evaluation 
of  the  low-frequency  ground  motion  response  characteris tics 
of  five  specific  prototype  MX  sites  and  were  shown  to  be  in 
good  agreement  with  the  available  finite  difference  calcula¬ 
tions.  In  particular,  the  scaling  laws  were  used  to  infer 
that  the  predominant  Rayleigh  wave  frequencies  at  these 
sites  will  depend  only  on  the  thickness  and  average  shear 
wave  velocity  of  the  surficial  soil  layer  and  will  be  in¬ 
dependent  of  yield  for  yields  greater  than  a  few  hundred 
kilotons . 

A  preliminary  analysis  of  the  Misers  Bluff  II-2  "800 
millisecond  anomaly"  was  presented  in  Section  2.5  where  it 
was  noted  that  both  the  spatial  dependence  of  the  anomaly  as 
well  as  the  inferred  axisymmetric  nature  of  its  source  are 
inconsistent  with  linear  superposition,  suggesting  that  the 
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observations  may  correspond  to  a  true  multiburst  effect. 
Furthermore,  it  was  shown  that  these  observations  can  be 
matched  remarkably  well  by  a  theoretical  model  which  con¬ 
siders  Rayleigh  waves  converging  on  the  center  of  the  source 
array  from  an  axisymmetric ,  impulsive  surface  source  acting 
at  a  radius  coincident  with  that  of  the  charge  locations. 

The  current  status  of  the  investigation  into  the 
valley  reverberation  problem  was  summarized  in  Section  III. 
This  summary  included  a  brief  overview  of  previous  experi¬ 
mental  and  theoretical  studies  as  well  as  a  detailed  de¬ 
scription  of  the  results  of  the  analysis  of  the  Misers  Bluff 
II  data.  These  latter  data  were  processed  using  a  variety 
of  frequency  and  direction-sensitive  filters  in  an  attempt 
to  isolate  the  source  locations  for  the  observed  late-time 
reverberation  arrivals.  Although  some  'Specific  arrivals  were 
identified  as  a  result  of  this  analysis,  the  interpretation 
of  the  data  was  complicated  due  to  the  limited  resolving 
power  of  the  seismic  arrays  and  the  complexity  of  the  valley 
structure.  However,  the  Misers  Bluff  II  data  were  found  to 
be  consistent  with  other  data  and  calculations  in  that  they 
indicate  that  the  late- time,  low-frequency  ground  motion  en¬ 
vironment  is  modified  by  the  presence  of  the  valley  bound¬ 
aries  . 

4.2  CONCLUSIONS 

The  analyses  summarized  above  support  the  following 
conclusions  concerning  the  characteristics  of  the  low-fre¬ 
quency  ground  motions  produced  by  near-surface  explosions. 

1.  The  seismic  source  function  for  airblast- 

induced  Rayleigh  waves  approaches  a  constant 
value  at  low  frequencies  (i.e.,  the  total  im¬ 
pulse)  and  decreases  as  frequency  squared 
-2 

(i.e.,  v  -jj  )  '  in  the  high-f recuency  limit. 


2.  At  low  frequencies  the  amplitude  of  the 
Rayleigh  wave  source  function  is  directly 
proportional  to  yield  (W)  and  independent 
of  the  site  geology.  At  high  frequencies 
the  spectral  amplitude  level  is  propor- 

j, 

tional  to  W3  and  increases  as  the  near¬ 
surface  shear  wave  velocity  increases. 

3.  For  the  prototype  MX  sites  considered  here, 
the  characteristic  corner  frequency  of  the 
Rayleigh  wave  site  response  function,  f  , 
is  given  approximately  by 


where  H  is  the  total  soil  thickness  at  the 
site  and  3  is  the  average  shear  wave  veloc¬ 
ity  of  the  soil  column. 

4.  For  yields  greater  than  a  few  hundred  kilo- 
tons  at  the  MX  prototype  sites,  the  predominant 
frequency  of  the  fundamental  mode  Rayleigh 
wave  motion  will  be  equal  to  f  ,  independent 

of  yield.  Moreover,  the  amplitude  of  that 
component  of  the  motion  will  scale 

it 

as  w3. 

5.  The  "800  millisecond  anomaly"  observed  in 
the  interior  of  the  MBII-2  source  array  has 
been  identified  as  a  converging  Rayleigh 
wave  produced  by  a  nearly  axisymmetric  load 
acting  on  a  small  surface  ring  centered  on 
the  charge  locations.  The  available  evidence 
suggests  that  it  is  a  true  multiburst  effect. 

With  regard  to  valley  reverberation  effects,  results 
studies  conducted  to  date  lead  to  the  following  conclusions. 
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The  available  data  confirm  that  the  pres¬ 
ence  of  the  valley  boundaries  can  signif¬ 
icantly  modify  the  late- time,  low-frequency 
ground  motion  environment  within  the 
valley.  In  particular,  the  duration 
of  this  component  of  the  motion  at  in¬ 
valley  stations  has  been  observed  to  be 
anomalously  long  by  factors  of  two  or  more 
with  regard  to  motion  measured  under  com¬ 
parable  conditions  outside  the  valley. 

The  characteristics  of  the  late-time  ground 
motion  data  recorded  from  the  MBII-1  and 
MBII-2  experiments  were  quite  different, 
despite  the  fact  that  they  were  detonated 
only  a  few  hundred  meters  apart  in  the  same 
valley.  This  indicates  that  the  valley  re¬ 
sponse  can  be  a  complex  function  of  the 
source  and  receiver  location  within  the 
valley  structure. 

Although  only  a  few,  highly  idealized  theo¬ 
retical  valley  calculations  have  been  per¬ 
formed  to  date,  they  have  consistently  pre¬ 
dicted  significant  boundary  effects  on  the 
late-time  ground  motion  environment  which 
are  in  qualitative  agreement  with  the  mea¬ 
sured  data. 
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